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A  MERCURY  MANOMETER  GRAVITY  GRADIOMETER 


by 

G .  Ian  Moore 
Frank  D.  Stacey 
Gary  J.  Tuck. 

Department  of  Physics 
University  of  Queensland 
Brisbane,  Australia 


ABSTRACT 

We  present  the  design  of  a  gravity  gradiometer  based  on  the  principles  of 
a  mercury  manometer.  The  gradiometer  consists  of  two  identical  manometers 
separated  vertically  by  1  meter.  The  pressure  difference  required  to  support 
the  mercury  columns  is  maintained  in  gas  chambers  above  pools  at  each  end  of 
the  mercury  columns  and  is  servoed  to  keep  the  lower  column  height  constant. 

A  change  in  the  vertical  gravity  gradient  leads  to  a  change  in  the  upper  column 
height.  The  column  heights  are  monitored  using  capacitance  micrometry.  Using 
the  lower  manometer  as  a  pressure  reference  and  servoing  to  a  constant  pressure 
removes  the  requirements  of  absolute  temperature  stability  and  dimensional 


stability  of  the  pressure  chambers. 
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Inrodurt ion : 

The  gravity  group  at  the  Physics  Department  of  the  University  of 
Queensland  is  currently  developing  a  gradiometer  based  on  the  principles 
of  a  mercury  manometer.  The  design  has  the  advantage  of  using  relatively 
simple  technology  and  has  a  target  sensitivity  useful  in  both  geophysical 
surveying  and  fundariiental  gravitational  experiments. 

Desinn  principle: 

The  principle  of  the  gradiometer  design  is  based  on  the  simple 
relationship  between  the  height  of  a  mercury  column  and  the  hydrostatic 
pressure  difference  between  the  two  mercury  surfaces.  By  way  of 
introduction,  consider  a  manometer  as  shown  in  Fig  1.  This  consists  of 
two  mercury  pools,  separated  by  0.5m  and  connected  by  a  fine  capillary 
tube.  The  pressure  difference  required  to  support  the  height  of  the 
mercury  column  is  maintained  in  the  gas  chambers  above  each  pool  and  is 
simply  given  by. 


^  ^  f’L  -  ^U  ^ 


where  p  is  the  density  of  the  mercury,  g  the  acceleration  due  to  gravity 


iiiul  il  the  column  height.  Generally  a  munti..ieter  is  used  to  delect  <h;iri)‘.'S 
in  the  pressure  difference,  howeiver  we  are  interested  in  keeping  the 
ptessun-  difference  constant  and  observing  small  changes  in  the  local 
gravitational  acceleration  by  way  of  a  change  in  the  column  height.  .Small 
changes  in  this  height  can  be  monitored  using  capaciteince  micrometry  to 
measure  the  changes  in  the  gap  between  each  mercury  surface  and  a 
stainless  steel  electrode  fixed  above  it  on  the  supporting  framework. 
Provided  the  requirement  of  constant  pressure  difference  can  be  met,  the 
maiiouieter  can  be  used  as  a  simple  vertical  axis  accelerometer  or 
gravimeter.  A  simple  calculation  shows  that  the  absolute  temperature  must 
hi-  inairit  aint:d  to  the  saioe  precision  as  that  desired  for  the  gravimeter. 

Kor  ail  instrument  capable  of  detecting  changes  as  small  as  10  of  g  (the 

c 

effect  due  to  a  gradient  of  1  E  over  a  distemce  of  1  metre),  this  is 
clearly  impractical.  Even  to  match  the  sensitivity  of  commonly  used 
gravimeters  (10  ^iGal)  requires  a  stability  of  order  micro  degrees  Kelvin 
which  is  extremely  difficult. 

The  requirements  for  a  very  stable  absolute  temperature,  a  high 
degree  of  dimensional  stability  in  the  gas  chambers  and  perfect  gas  seals 
make  a  useable  gravimeter  based  on  this  principle  impractical.  Fortunately 
these  problems  can  be  greatly  reduced  in  the  design  of  a  gradiometer. 

The  gradiometer  is  formed  by  coupling  two  single  manometers  as  shown 
in  Fig  "g.  The  lower  chambers  of  each  manometer  are  connected  by  a  small 
lube  as  are  the  upper  chambers.  The  pressure  difference  across  the  pairs 
of  chambers  is  servoed  by  means  of  an  adjustable  bellows  to  maintain  the 
lower  column  height  constant.  Any  change  in  the  gravity  difference  acting 
on  the  two  columns  due  to  a  change  in  vertical  gravity  gradient  will  lead 
to  a  small  change  in  the  height  of  the  upper  column. 


Since  the  lower  manometer  now  forms  a  precise  pressure  refertuice  tlic 
absolute  temperature  of  the  instrument  is  no  longer  important.  It  is  only 
necessary  that  temperature  gradients  along  the  length  of  the  instrument  be 
kept  small.  Also,  since  the  pressure  chambers  of  the  manometers  are 
directly  coupled,  any  slight  leak  or  dimensional  change  in  any  gas  chamber 
becomes  a  common  mode  signal  and  is  removed  by  the  pressure  servo. 

Design  details: 

I .  .Sensitivity: 

As  a  basis  for  the  design  we  take  a  target  sensitivity  of  1  F.,  equal 
to  that  generally  accepted  as  useful  in  terms  of  geophysical  surveying. 

For  a  column  height  of  0.5  m  and  a  separation  of  1  m  between 
corresponding  pools  of  the  coupled  manometers  as  shown  in  Fig  2,  a  change 
in  gradient  of  1  E  gives  rise  to  a  change  in  the  upper  coluiun  height  of 
0.05  nm  or  half  of  this  as  the  detectable  change  in  the  level  of  a  pool. 
This  displacement  comes  close  to  the  observable  limit  of  the  capacitance 
micrometry  technique  used.  This  technique  will  be  discussed  in  more  detail 
in  part  IV  of  this  section. 

The  mercury  pools  act  as  pistons  as  the  column  heights  change, 
varying  the  gas  volume  and  hence  the  supporting  pressure.  This  effect 
tends  to  stabilize  the  column  heights  against  perturbations  and  therefore 
would  reduce  the  sensivity  of  a  single  manometer  to  gravitational  changes. 
However,  in  a  differential  instrument  this  problem  does  not  arise  because 


3. 


we  are  concerned  with  height  differences  of  mercury  columns  subjected  to  a 
(ronunon  pressure  difference. 

I I .  Thermal  stabilization: 

There  are  two  requirements  of  the  thermal  stabilization  of  the 
instrument.  First,  to  maintain  a  stable  absolute  temperature  and  second, 
to  maintain  a  constant  (preferably  zero)  temperature  gradient  along  the 
lengtli  of  the  instrument. 

The  thermal  expansion  of  the  mercury  dominates  the  stability  problem 

because  we  are  concerned  with  thermal  changes  in  its  density  relative  to 

the  thermal  expansion  of  the  stainless  steel  framework  of  the  instrument. 

-4  -1 

The  relevant  parameter  is  (a,,  -  a  )  where  s*  2  x  10  K  is  the 

Hg  ss  Hg 

-5  -1 

volume  expansion  coefficient  for  mercury  and  a  a  1.7  x  10  K  is  the 

SS 

linear  coefficient  for  the  stainless  steel.  If  it  were  necessary  to 
obtain  absolute  accuracy  of  1  part  in  10^^  in  a  single  manometer  then 
absolute  temperature  stability  of  0.5  fjK  would  be  required  emd  this  cannot 
be  realized.  In  a  differential  mode  we  formally  require  this  precision  in 
the  temperature  difference  between  manometers  and  the  thermal 
stabilization  system  described  below  is  targetted  on  this  accuracy.  This 
looks  a  difficult  target  but  we  have  an  internal  check  on  strong 
tt.-iuperuture  gradients,  because  levels  in  all  four  mercury  pools  are 
monitored.  The  instrument  can  actually  measure  the  thermal  expansion  of 
the  mercury  (relative  to  the  plastic  pools  and  the  stainless  steel)  in 
each  manometer  independently,  although  it  is  not  clear  that  this  will  be 
directly  useful  and  further  developments  to  reduce  thermal  sensitivity  are 
uniior  considerat  ion. 
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There  is  a  reisidual  thermal  problem  if  the  two  manuiuet  i^rr.  are 
slightly  different  dimensionally,  but  assuming  matched  lengths  of  the 
stainless  steel  spacing  rods  to  better  than  1  part  in  10^,  which  is  <,*asily 
achieved,  O.Ofi  K  absolute  temperature  stability  suffices  to  avoid  this 
problem. 

The  temperature  stabilization  system  which  we  propose  to  use  to 

achieve  the  above  requirements  is  shown  in  Fig  3.  This  consists  of  a 

double  circulated  water  jacket,  the  temperature  of  the  water  being 

controlled  at  a  point  just  prior  to  the  entry  into  the  instrument  jacket 

by  means  of  a  heater  and  thermistor  feedback.  The  water  first  circulates 

down  the  inner  jacket  and  then  back  up  the  outer  jacket.  The  outer  casing 

is  covered  with  closed  ciill  insulation  and  the  wall  between  the  inner  and 

outer  jackets  is  also  well  insulated.  The  instrument  itself  is  housed 

inside  an  inner  casing.  In  principle  it  should  be  possible  to  achieve 

_6 

temperature  gradients  of  less  than  2  x  10  K  over  the  length  of  the 

instrument  at  the  inner  water  jacket  for  an  external  temperature 

controlled  to  within  1  K.  The  final  passive  shield  consisting  of  a  layer 

of  insulation  over  the  copper  instrument  case  should  reduce  this  to  less 

-7 

than  the  required  5  x  10  K  over  the  length  of  the  instrument.  A 
completely  passive  thermal  shield  (relying  on  thermal  conduction  rather 
than  a  circulating  media)  which  would  achieve  a  similar  result  is  totally 
impractical  because  of  the  large  surface  area  to  cross  section  ratio  of 
the  shielding  tubes. 

III.  .Servo  Bellows: 

The  servo  bellows  is  currently  driven  in  two  stages.  Firstly  with  a 
coarse  motor  driven  differential  thread  and  secondly  with  a  piezo-stack 
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for  fine  control.  This  is  a  temporary  arranj^eraent  and  the  final  .servo  is 
to  be  driven  using'  an  "inch-worm”  device  winch  will  give  fine  control  ovi.-r 
a  rang!'  of  6iiuii  and  remove  the  problems  of  backlash  etc  associati-d  with  tin- 
mechanical  drive.  The  large  force  exerted  by  the  bellows  on  its  drive  is 
overcome  by  plairing  the  servo  bellows  in  a  chamber  pressurised  to  very 
nearly  the  same  pressure  as  the  lower  manometer  chambers  to  wliich  the 
bellows  is  connected. 


DcttM't  ion  System: 


Between  each  pool  and  its  corresponding  gas  chamber  is  an  arrangement 

of  stainless  steel  electrodes,  detailed  in  Fig  4.  This  forms  a  fixed 

capacitance  gap  against  which  the  capacitance  between  the  mercury  surface 

and  tile  central  electrode  is  compared.  The  technique  of  capacitance 

micrometry^  is  described  with  reference  to  Fig  5.  The  upper,  fixed 

electrode  and  the  mercury  are  excited  in  antiphase  with  3  kHz  (3  Vpeak) 

sigiuils  derived  from  a  switchable  ratio  transformer.  The  resulting  3  kHz 

signal  on  the  central  electrode  is  detected  synchronously  with  the 

excitation  signal.  The  ratio  transformer  switch  setting  is  adjusted  until 

the  detected  signal  is  a  null.  The  switch  setting  then  gives  a  direct 

reading  of  the  ratio  of  the  two  gaps.  The  mercury  gap  is  then  calculated 

from  this  ratio  and  the  known  fixed  gap.  Only  the  first  five  digits  of 

the  ratio  are  obtained  from  the  ratio  transformer;  three  more  digits  are 

obtaineil  by  measuring  the  out  of  balance  signal  at  the  final  ratio 

I  rails  former  setting.  Thus  with  capacitance  gups  of  0.2  mm  the  detector 
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sensitivity  is  better  than  10  of  this  or  <0.02  nm.  This  corresponds 

approximately  with  the  expected  change  in  a  single  gap  for  a  change  in 


grail ient  of  I  K. 


An  automatic  ratio  transformer  bridge  samples  the  four  channc-ls 
corresi<onding  to  the  upper  ;uul  lower  pools  of  each  manoiueter  und<i' 
computer  control.  The  computer  controls  the  pressure  servo  feedba<’k  loop 
and  provides  automatic  readout  of  the  gaps  and  column  heights. 

V .  Meroui V  Pool  Floats: 

Potentially  the  most  significant  problem  and  one  which  is  currently 
being  addressed  is  that  of  rippling  of  the  mercury  surfree  when  the 
instrument  is  vibrated.  This  rippling  causes  two  problems:  First,  slight 
impurities  in  the  mercury  make  the  mercury  stick  to  the  stainless  steel 
electrodes  when  the  gap  is  very  small.  Second,  and  more  importantly  from 
a  fundamental  view  point  is  that  this  rippling  causes  a  bias  in  the 
measured  capacitance  gaps.  This  is  due  to  the  fact  that  the  average 
reciprocal  gap  is  measured  and  this  is  biased  from  the  reciprocal  of  the 
av«;rage  gap  when  the  mercury  sut  face  is  not  flat  as  shown  Fig  6. 

To  overcome  this  problem,  we  intend  to  use  stainless  steel  floats  on 
the  mercury  pools.  These  will  be  constrained  by  flat  stainless  steel 
springs  to  prevent  them  drifting  to  one  side.  The  proposed  arrangement  is 
shown  in  Fig  7.  The  springs  are  etched  from  O.OOi"  sheet  (Fig  8)  and  ari* 
designed  to  have  a  large  compliance  to  motions  perpendicular  to  their 
plane  and  retain  a  high  stiffness  to  motions  in  their  plane.  The  floats 
are  rebated  in  order  to  avoid  possible  interference  between  the  spring 
leaves  and  their  lower  faces. 

The  floats  rest  on  supports  until  the  mercury  pools  are  filled  and 
the  spring  anchors  are  machined  so  that  the  desired  working  gaps  are 
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obtiiiiifd  with  zero  extension  of  the  springs.  This  is  important  for  two 
leasons.  First,  the  stiffness  of  the  springs  rises  rapiclly  with  large 
extensions  and  very  compliant  springs  are  required  if  the  sensitivity  of 
the  last  ruint.*nt  is  not  to  be  reduced.  Secon.l,  the  natural  bouyancy  level  of 
the  floats  must  not  be  significantly  altered  by  the  springs.  In  the 
absence  of  springs  the  floats  would  always  float  at  the  same  level  in  the 
mercury  regardless  of  the  gravitational  acceleration.  Since  the  springs 
may  exert  a  no:  -zero  force  in  the  vertical  direction,  this  will  no  longer 
be  true  as  illustiated  in  Fig  U. 


where. . . 


M 


ks 


--(2) 


6  is  the  flotation  level  of  a  float 
M  is  the  mass  of  a  float 
A  is  the  cross  sectional  area  of  a  float 

p,,  and  p  are  the  desities  of  mercury  and  stainless  steel 
'^Hg  '^ss 

s  is  the  effective  operating  extension  of  the  paired  springs 
k  is  the  effective  spring  constant  of  the  paired  springs 


Hence  if  the  gravitational  field  changes  (for  example  with  a  gain  in 

altituile)  the  level  at  which  the  floats  ride  in  the  mercury  will  change 

for  each  individual  float,  thereby  voiding  the  measurements.  Assuming  a 

4 

change  in  gravity  of  one  part  in  10  and  requiring  the  natural  flotation 
level  of  the  floats  to  remain  constant  within  the  tolerance  implied  by  the 
target  sensivity  gives  the  following  condition, 


ks 


X  10 


<  0.2  X  10  m 


- (3) 
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Assuming  the  operating  extension  of  the  spring  can  be  kept  to  less  than 
0. 1  mm  this  gives  an  upper  limit  on  the  spring  constant  of  0.14  Nm  So 
the  springs  must  be  extremely  compliant,  at  least  over  the  very  limited 
range  of  extensions  expected  during  operation. 

References: 

1:  STACEY,  F.D.,  HYNN,  J.M.W.,  LITTLE,  E.C.  and  CROSKELL,  C. 

Displacement  and  tilt  transducers  of  140  dB  range.  J.  Sci . 
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Fiirun*  captions: 


Figure  1. 

Basil’  schematic  of  a  single  mercury  manometer. 

Figure  2. 

Schematic  of  a  double  manometer  showing  interconnections  of  pressure 
chambers  and  servo  bellows. 

Figure  2. 

Proposed  temperature  shield  consisting  of  a  double,  insulated,  circulated 
water  Jacket  with  heater  and  thermistor  feedback.  The  instrument  is 
housed  inside  a  final  passive  thermal  shield  consisting  of  a  layer  of 
insulation  around  a  heavy  copper  casing.  The  whole  assembly  operates  in 
an  environment  stabilized  to  within  1  K. 

Figure  4. 

Stainless  steel  electrode  configuration  used  to  meetsure  small  changes  in 
the  height  of  the  mercury  column  by  the  technique  of  capacitance 
micrometry. 


Figure  5. 

General  schematic  of  electronics  used  for  capacitance  micrometry. 


Figure  6. 

Oiagramatic  view  of  rippling  on  the  mercury  surface.  In  general  the 
average  reciprocal  gap  measured  by  the  capacitance  micrometry  is  less  than 
the  reciprocal  of  the  average  gap. 
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Figure  7. 

Proposed  arrangements  of  stainless  steel  boats  and  constraining  springs 
intended  to  overcome  the  problem  of  rippling  of  the  mercury  surface. 

Figure  8. 

Detail  of  spring  shape.  These  are  etched  from  0.001"  stainless  steel  sheet 
and  have  a  very  high  compliance  in  a  direction  perpendicular  to  the  plane 
of  the  spring. 

figure  9. 

The  effect  of  a  spring  on  the  natural  bouyancy  of  a  float.  The  flotation 
level  is  no  longer  independent  of  g  and  the  spring  constant  must  be  very 


small. 
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TITLE  OF  PAPER:  GRADIO  Project:  High  Sensitivity  Electrostatic 
For  Spaceborne  Gradiometry 

SPEAKER:  Georges  Balmino 
QLESTIONS  AND  CQ^LMEN•TS: 

1.  Question:  Ho  Jung  Paik 

Could  you  go  over  how  you  get  out  of  the  cross  coraponents? 
You  have  T'.t? ,  but  do  you  know  Af? 

Response : 

We  use  star  trackers  and  gyros  for  that. 

2.  Question :  Jean-Paul  Richard 

Accuracy  of  star  trackers? 

Response : 


Accelerometer 


Star  trackers  accuracy  =  0.01  sec  of  arc/sec 
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CUBE  MOTION  DETECTION 


Translation  :  S  -  5  ‘  =  4  ax 


ACCELEROMETER  RESOLUTION 


4 


000 


digitalization 

alignment  of  the  accelerometers 
scale  factor  equalization 


I 


) 


GO 

M 

C0 

a; 

0) 

a 

N  2 

w  S 

•  o 
o  o 

d 


pCt 

o 

d 

d 

o 


CV2 

I 


d 

o 


00  q; 

a  § 


•  « 


CM 

I 


CM 

I 

00 

a 

•  • 

CM 

00  ^ 

I 

O  _ 

^  o 

"•H 

^  fl' 

>5: 

fl;  > 

d  < 


^  d 


d  GO 

^  i3 

o  TJ 

jl 

o  P, 

*te 


a; 


0? 

jL4 

d 

GO 

d 

0) 

a 

0) 

d 

o 


o 

00  < 
0)  ^ 
K  > 


d 


(  10  “J  Jii.u.  for  a  1  m  base  line  ) 

necessary  data  for  calibration  processing 

(  simple  dedicated  algorithms  ) 


ACCELEROMETER  SCALE  FACTOR  AND  CROSS-COUPLING  COEFFICIENTS 
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ACCELEROMETERS  :  SENSITIVITY  ERRORS 
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SENSITIVITIES  TO  DISPLACEMENTS 


in  a  spherical  approximation  and  local 
orbital  axes  : 


^  Gravity  gradient  components  are  not 
sensitive  to  horizontal  displacements 


TRANSFORMATION  OF  COORDINATES  : 
SPACECRAFT  TO  LOCAL  ORBITAL  AXES 


Uncertainty  due  to  errors  dSj  ,  50l  • 
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^^>Earth  pointing  (  0j,0l,0r«1  )  .* 

T,j  are  not  sensitive  to  60l 

(rotation  about  the  vertical  axis) 

OAltitude  200  Km  : 
b0j=  1  arc  sec 
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and  T22  are  useful  for  navigation  (  Estimation  of  6r  ) 


GRADIOMETRY  THROUGH  DIFFERENTIAL  ACCELEROMEI HY 


GRAVITY  GRADIENT  ,  ORIENTATION  ,  POSITION 


Fifteenth  Gravity  Gradiometer  Conference 
United  States  Air  Force  Academy 
Colorado  Springs,  Colorado 

CONFERENCE  AGENDA 


Tuesday,  10  February  1987 

1900  -  2200  -  Pre-Conference  Get-Together  at  Hilton  Inn 

Early  Registration 

Wednesday,  11  February  1987 

0700  -  Depart  Hilton  Inn  for  Fairchild  Hall 

0730  -  Registration  -  3rd  floor  Fairchild  Hall,  South  End 

0745  -  Welcome/Introduction  -  Capt  Terry  J-  Fundak 

0815  -  Presentation  by  Dr.  Georges  Balmino  of  the  ONERA  (Office  National 
d’Etudes  et  de  Recherches  Aerospatiales) . 

"GRADIO  Project:  A  SGG  Mission  Based  on  Microaccelerometers” 

0845  -  Presentation  by  Dr.  G.  Ian  Moore  of  the  University  of  Queensland. 

"A  Mercury  Manometer  Gravity  Gradiometer" 

0900  -  Presentation  by  Mr.  Ernest  H.  Metzger  of  Bell  Aerospace  Textron. 

"Bell  Aerospace  Gravity  Gradiometer  Survey  System  (GGSS)  -  Program 
Review" 

0925  -  Presentation  by  Dr.  Frank  J.  van  Kann  of  the  University  of  Western 
Australia. 

"A  Prototype  Superconducting  Gravity  Gradiometer  for  Geophysical 
Exploration" 

0952  -  Presentation  by  Dr.  Warren  G.  Heller  of  The  Analytic  Sciences  Corp. 

"Gravity  Gradiometer  Survey  System  (GGSS)  Data  Processing  and  Data 
Use" 

1016  -  Break 

1035  -  Presentation  by  Mr.  A1  Jlrcitano  of  Bell  Aerospace  Textron. 


"Self-Gradient  Calibration  of  the  GGSS  in  a  C-130  Aircraft" 


1058  -  Presentation  by  Dr.  Sam  C.  Bose  of  Applied  Sciences  Analytics,  Inc. 

"Gravity  Gradiometer  Data  Processing  Using  the  Karhunen-Loeve 
Method" 

1120  -  Presentation  by  Mr.  David  M.  Gleason  of  the  Air  Force  Geophysics 
Laboratory. 

"Numerically  Deriving  the  Kernels  of  an  Integral  Predictor  Yielding 
Surface  Gravity  Disturbance  Components  from  Airborne  Gradient  Data" 

1130  -  Presentation  by  Mr.  A1  Jircitano  of  Bell  Aerospace  Textron. 

"Stage  II  Simulation  Results  Using  the  NSWC  Synthetic  Gravity 
Field" 

1150  -  Depart  Fairchild  Hall  for  USAFA  Noncommissioned  Officers' 

(NCO)  Club 

1200  -  Lunch  -  USAFA  NCO  Club 

1245  -  Depart  NCO  Club  for  Fairchild  Hall 

1330  -  Presentation  by  Dr.  Richard  H.  Rapp  of  Ohio  State  University. 

"Gradient  Information  in  New  High  Degree  Spherical  Harmonic 
Expansions" 

1354  -  Presentation  by  Mr.  John  J.  Graham  of  the  Defense  Mapping  Agency 
Aerospace  Center. 

"The  Effect  of  Topography  on  Airborne  Gravity  Gradiometer  Data" 

1357  -  Presentation  by  Mr.  Mike  Sideris  of  the  University  of  Calgary. 

"Effect  of  Terrain  Representation,  Grid  Spacing,  and  Flight  Altitude 
on  Topographic  Corrections  for  Airborne  Gradiometry" 

1417  -  Presentation  by  Dr.  Rene  Forsberg  of  Geodetic  Institute  (Denmark) 
(Currently  at  the  University  of  Calgary,  Canada). 

"Topographic  Effects  in  Airborne  Gravity  Gradiometry" 

1434  -  Presentation  by  Dr.  Alan  H.  Zorn  of  Dynamics  Research  Corporation. 

"Observability  of  Laplace's  Equation  Using  a  Torsion-Type 
Gravity  Gradiometer" 

1510  -  Break 

1530  -  Presentation  by  Dr.  Carl  Bowln  of  Woods  Hole  Oceanographic  Institute 

"Ratios  of  Gravity  Gradient,  Gravity,  and  Geold  for  Determination  of 
Crustal  Structure" 


1550  -  Presentation  by  Dr.  Rene  Forsberg  of  Geodetic  Institute  (Denmark). 

"Combining  Gravity  Gradiometry  with  Other  Exploration 
Methods  for  Geophysical  Prospecting" 

1600  -  Presentation  by  Dr.  Rene  Forsberg  of  Geodetic  Institute  (Denmark). 

"Computation  of  the  Gravity  Vector  from  Torsion  Balance  Data  in  S. 
Ohio" 

1615  -  Presentation  by  Dr.  Hans  Baussus  vcn  Luetzow  of  the  U.S.  Army 
Engineer  Topographic  Laboratories. 

"Estimation  of  Gravity  Vector  Components  from  Bell  Gravity  Gradiometer 
and  Auxiliary  Data  under  Consideration  of  Topography  and  Associated 
Analytical  Upward  Continuation  Aspects" 

1635  -  Depart  Fairchild  Hall  for  the  Hilton  Inn 

1700  -  Reception  -  Hilton  Inn 

Thursday,  12  February  1987 

0700  -  Depart  Hilton  Inn  for  Fairchild  Hall 

0755  -  Presentation  by  Dr.  M.  Vol  Moody  of  the  University  of  Maryland. 

"Development  of  A  Three-Axis  Superconducting  Gravity  Gradiometer 
and  a  Six-Axis  Superconducting  Accelerometer" 

0835  -  Presentation  by  Dr.  Bahram  Mashhoon  of  the  University  of  Missourl- 
Columbla. 

"The  Gravitational  Magnetic  Field  of  the  Earth  and  the  Possibility 
of  Measuring  It  Using  an  Orbiting  Gravity  Gradiometer" 

0905  -  Presentation  by  Dr.  Ho  Jung  Paik  of  the  University  of  Maryland. 

"Tests  of  General  Relativity  in  Earth  Orbit  Using  a  Superconducting 
Gravity  Gradiometer" 

0928  -  Presentation  by  Dr.  Dave  Sonnabend  of  Jet  Propulsion  Laboratory. 

"Magnetic  Isolation  -  Closing  the  Loop" 

0941  -  Presentation  by  Dr.  Dan  Long  of  Eastern  Washington  University. 
"Laboratory  G(R)  Experiment  -  Progress  Report" 


1004 


Break 


1030  -  Cheyenne  Mountain  Complex  Overview 
Briefing  by  Maj  Bill  Carver,  L’SAF 
(Chief,  NORAD  Presentations  Division). 

1110  -  Form  Groups  A  &  B 

1115  -  Depart  Fairchild  Hall  for  USAFA  NCO  Club 
1130  -  Lunch  -  USAFA  NCO  Club 

1200  -  Depart  USAF  Academy  for  Falcon  Air  Force  Station 

1245  -  Arrive  Falcon  AFS  for  briefing  on  2nd  Space  Wing 

Tour  of  the  Consolidated  Space  Operations  Center  (CSOC) 

(Group  A) 

1415  -  Depart  CSOC  for  Cheyenne  Mountain  Complex  (CMC) 

1500  -  Arrive  CMC 

1505  -  Security  in-processing  and  process  through  metal  detector 
1525  -  Travel 

1530  -  Tour  NORAD  Command  Post 
Tour  Industrial  Area 

1620  -  Travel/question  and  answer  session 
1630  -  Depart  for  Hilton  Inn 

1715  -  Arrive  Hilton  Inn 

(Group  B) 

1415  -  Depart  USAF  Academy  for  Peterson  Air  Force  Base  (AFB) 

1445  -  Arrive  Peterson  AFB  museum 

1600  -  Depart  Peterson  AFB  for  Hilton  Inn 

1630  -  Arrive  Hilton  Inn 


Friday,  13  February  1987 

0800  -  Tour  of  JILA,  Boulder,  Colorado 


Papers  Included  in  VOLUME  I  of  the  Conference  Proceedings 


1.  *Dr.  Georges  Balmlno,  C.N.E.S. /Bureau  Gravlmetrlque  International,  France 
Dr.  Alain  Bernard,  ONERA  (Office  National  d'Etudes  et  de  Recherches 

Aerospat iales ,  Rrance) 

Dr.  Pierre  Touboul,  ONERA,  France 

"GRADIO  Project  :  A  SGG  Mission  Based  on 
Microaccelerometers" 


2.  *Dr.  G.  Ian  Moore,  University  of  Queensland,  Australia 
Dr.  Frank  D.  Stacey,  University  of  Queensland,  Australia 
Dr.  Gary  J.  Tuck,  University  of  Queensland,  Australia 
Dr.  Barry  D.  Goodwin,  University  of  Queensland,  Australia 

"A  Mercury  Manometer  Gravity  Gradiometer" 


3.  Mr.  Louis  L.  Pfohl,  Bell  Aerospace  Textron 
*Mr.  Ernest  Metzger,  Bell  Aerospace  Textron 

"Bell  Aerospace  Gravity  Gradiometer  Survey 
System  (GGSS)  -  Program  Review” 


4.  *Dr.  Frank  J.  van  Kann,  et  al.  University  of  Western  Australia 

"A  Prototype  Superconducting  Gravity 
Gradiometer  for  Geophysical  Exploration” 


5.  *Dr.  Warren  G.  Heller,  The  Analytic  Sciences  Corporation 

"Gravity  Gradiometer  Survey  System  (GGSS) 
Data  Processing  and  Data  Use” 


6.  Dr.  W.  John  Hutcheson,  Bell  Aerospace  Textron 

(Paper  presented  by  Mr.  Al  Jlrcitano  of  Bell  Aerospace  Textron) 

"Self-Gradient  Chlibration  of  the  GGSS 
in  a  C-130  Aircraft” 


7.  *Dr.  Sam  C.  Bose,  Applied  Science  Analytics,  Inc 

Mr.  Glenn  E.  Thobe,  Applied  Science  Analytics,  Inc 

"Gravity  Gradiometer  Data  Processing  Using 
the  Karhunen-Loeve  Method" 


*  Denotes  Speaker  at  Conference 


8*  *Mr .  David  M.  Gleason,  Air  Force  Geophysics  Laboratory 


"Numerically  Deriving  the  Kernels  of  an  Integral 
Predictor  Yielding  Surface  Gravity  Disturbance 
Components  from  Airborne  Gradient  Data" 


9.  Dr.  W.  John  Hutcheson,  Bell  Aerospace  Textron 

(Paper  presented  by  Mr.  A1  Jircltano  of  Bell  Aerospace  Textron) 

"Stage  II  Simulation  Results  Using 
the  NSWC  Synthetic  Gravity  Field” 


10.  *Dr.  Richard  H.  Rapp,  Ohio  State  University 

"Gradient  Information  in  New  High 
Degree  Spherical  Harmonic  Expansions 


*  Denotes  Speaker  at  Obnference 


Papers  Included  In  VOLUME  IL  of  the  Conference  Proceedings 


1.  *Mr.  John  J.  Graham,  Defense  Mapping  Agency  Aerospace  Cbnter 

Mr.  Joseph  L-  Toohey,  Defense  Mapping  Agency  Aerospace  Center 

"The  Effect  of  Topography  on  Airborne 
Gravity  Gradiometer  Data" 


2.  Dr.  Klaus-Peter  Schwarz,  University  of  Chlgary,  Cknada 
*Mr.  M.G.  Siderls,  University  of  Calgary,  Canada 
Dr.  I.N.  Tzlavos,  University  of  Qilgary,  Chnada 
(Dr.  Tziavos  on  leave  from  the  University  of  Thessaloniki,  Greece) 

"Effect  of  Terrain  Representation,  Grid  Spacing,  and 
flight  Altitude  on  Topographic  Corrections  for 
Airborne  Gradlometry" 


3.  *Dr.  Rene  Forsberg,  Geodaetlsk  Institut,  Denmark 

"Topographic  Effects  in  Airborne  Gravity  Gradlometry" 


4.  *Dr.  Alan  H.  Zorn,  Dynamics  Research  Garporation 

"Observability  of  Laplace’s  Equation  Using 
a  Torsion-Type  Gravity  Gradiometer" 


5.  *Dr.  Carl  Bowin,  Woods  Hole  Oceanographic  Institute 

"Ratios  of  Gravity  Gradient,  Gravity,  and  Geoid 
for  Determination  of  Crustal  Structure" 


6.  Dr.  Anthony  A.  Vassiliou,  University  of  Calgary,  Canada 

(Paper  presented  by  Dr.  Rene  Forsberg,  Geodaetisk  Institut,  Denmark) 

"Combining  Gravity  Gradiometry  with  other 
Exploration  Methods  for  Geophysical  Prospecting 


7.  Dr.  D.  Arabelos,  University  of  Thessaloniki,  Greece 

Mr.  Christian  Tscherning,  Geodaetisk  Institut,  Denmark 

(Paper  presented  by  Dr.  Rene  Forsberg,  Geodaetisk  Institut,  Denmark) 

"Computation  of  the  Gravity  Vector  from  Torsion 
Balance  Data  in  Southern  Ohio” 


*  Denotes  Speaker  at  Conference 


8-  *Dr.  Hans  Baussus  von  Luetzow,  US  Army  Engineer  Topographic  Laboratory 

"Estimation  of  Gravity  Vector  Components  from  Bell  Gravity 
Gradiometer  and  Auxiliary  Data  under  Consideration  of 
Topography  and  Associated  Analytical  Upward  Continuation 
Aspects" 


9.  Dr.  H.  A.  Chan,  University  of  Maryland 
Dr.  Q.  Kong,  University  of  Maryland 
*Dr.  M.  Vol  Moody,  University  of  Maryland 
Dr.  H.  J.  Paik,  University  of  Maryland 
Mr.  J.  W.  Parke,  University  of  Maryland 

"Development  of  a  Three-Axis  Superconducting  Gravity 
Gradiometer  and  a  Six-Axis  Superconducting  Accelerometer" 


10.  *Dr.  Bahram  Mashhoon,  University  of  Missouri- Columbia 

"The  Gravitational  Magnetic  Field  of  the  Earth  and 
the  Possibility  of  Measuring  it  Using  an  Orbiting 
Gravity  Gradiometer" 


11.  *Dr.  Ho  Jung  Paik,  University  of  Maryland 

"Tests  of  General  Relativity  in  Earth  Orbit 
Using  a  Superconducting  Gravity  Gradiometer" 


12.  *Dr.  Dave  Sonnabend,  Jet  Propulsion  Laboratory 

Mr.  A.  Miguel  San  Martin,  Jet  Propulsion  Laboratory 

"Magnetic  Isolation- Closing  the  Loop" 


13.  *Dr.  Dan  Long,  Eastern  Washington  University 

"Laboratory  G(R)  Experiment  -  Progress  Report" 


*  Denotes  Speaker  at  Conference 


GRADIO  PROJECT: 

A  SGG  MISSION  BASED  ON  MICROACCELEROMETERS 


by 

Dr.  Georges  Balmino 

C. N. E. S . /Bureau  Gravlmetrique  International 
18  Ave  Edouard  Berlin 
31055  Toulouse  Cedex 
FRANCE 

Dr.  Alain  Bernard 
Dr.  Pierre  Touboul 

Office  National  d' Etudes  et  de  Recherches  Aerospatiales 

BP  72 

92322  Chatillon  Cedex 
FRANCE 


ABSTRACT 

The  status  of  the  satellite  gravity  gradiometry  project  is  reviewed.  Since 
the  first  ideas  in  1980,  technological  solutions  have  ripened  and  a  configura¬ 
tion  composed  of  eight  cubic  electrostatic  microaccelerometers  is  proposed  which 
should  guarantee  a  signal  detection  limit  of  10~2  to  10“^  Eotvos.  Two  basic 
systems  are  proposed  to  fly  the  instrument  :  one  is  a  dedicated  satellite  on 
which  a  permanent  calibrating  device,  actually  part  of  the  gradiometer,  would 
be  implemented;  the  other  would  consist  of  flying  the  instrument  in  one  of  the 
NASA  projects,  the  GRM  drag-free  spacecraft,  where  it  would  be  suspended  in  the 
double  stage  DISODS  system.  A  laboratory  model  of  the  cubic  accelerometers  is 
also  presented. 


GRADIO 


A  SGG  MISSION  BASED  ON  MICROACCELEROMETERS 


BERNARD  A.  (0NERA>  ChatUlon  s/Bagneuxi  Prance) 
TOUBOUL  P.  (ONERA.  id.) 

BALMINO  G.  <CNES.  Toulouse>  France) 


15th  GRAVITY  GRAOIOMETRY  CONFERENCE 
Feb.  11-13.  1987 
(Colorado  Springs  -  USA) 


RESOLUTION  -  100  KM) 


must  be  controlled  and  estimated 


GRADIO  ACCELEROMETER  ;  SPECIFICATIOHS 


O  low  coupling 
O  occic  fcctorc  mctehing 


TRANSFORMATION  FROM  SPATIAL  DOMAIN  TO  KARHUNEN-LOEVE  DOMAIN 


DEFINITION  OF  TRANSFORMATION  MATRICES 


APPLIED  SCIENCE  ANALYTICS.  INC. 


ABOUT  THE  GRAVITY  GRADIOMETRY  CONFERENCE 


The  First  Gravity  Gradiometry  Conference  was  held  at  the  Air  Force  Cambridge 
Research  Laboratory  (AFCRL,  now  AFGL)  in  1973.  Its  purpose  was  to  provide 
a  forum  to  evaluate  and  compare  the  efforts  of  three  vendors  (Charles  Stark 
Draper  Lab,  Hughes  Research  Lab  and  Bell  Aerospace  Textron)  in  still-emerging 
areas  of  gravity  gradiometry.  About  15  people  attended,  most  of  them  from  the 
companies  mentioned  above  or  the  Terrestrial  Sciences  Division  at  AFCRL.  In 
contrast,  the  1987  Conference  had  a  guest  list  of  70  plus  attendees,  with 
participation  from  academia  (foreign  and  domestic),  private  industry  and 
government.  The  papers  presented  were  not  restricted  to  gradiometry  alone. 
Indeed,  the  scope  of  this  annual  event  has  broadened  considerably  since  1973. 

With  the  exception  of  the  first  two  conferences,  all  the  others  have  been 
held  at  the  US  Air  Force  Academy  in  Colorado  Springs,  Colorado.  The  Geodesy 
and  Gravity  Branch  of  the  Earth  Sciences  Division  of  the  Air  Force  Geophysics 
Laboratory  (AFGL),  Hanscom  AFB,  Massachusetts,  has  always  organized  the  event, 
which  usually  takes  place  around  the  second  week  in  February.  This  trend  is 
expected  to  continue. 

If  you  are  not  already  on  our  mailing  list  and  would  like  to  attend  the 
1988  Conference,  or  if  you  have  any  questions,  please  write  to: 

Ms  Claire  McCartney 
AFGL/LW 

Hanscom  AFB,  MA  01731 

Due  to  space  constraints,  we  restrict  the  size  of  our  Conferences  to  about  75 
people.  Attendance  will  generally  be  on  a  "first-come,  first-served  basis"  once 
the  completed  registration  forms  are  returned  to  us.  We  shall  mail  these  forms 
later  this  year. 

While  we  have  a  limited  number  of  copies  of  the  proceedings  for  non-attendees 
of  the  1987  Conference,  copies  of  proceedings  for  prior  years  are  not  available. 
Also,  we  appreciate  any  comments  or  suggestions  you  may  have  regarding  this 
document . 


ABOUT  THESE  PROCEEDINGS 


Due  to  the  large  number  of  papers  presented  at  the  Conference,  I  have 
divided  the  proceedings  into  two  manageable  volumes.  At  the  beginning  of 
each  volume  is  a  list  of  all  papers  contained  in  both  volumes,  in  actual 
order  of  presentation  at  the  Conference.  This  is  also  the  sequence  of  the 
published  papers  within  these  proceedings. 

For  the  sake  of  completeness,  both  volumes  contain  the  Attendees  List, 
Conference  Agenda,  Lists  of  Papers,  Conference  History,  Acknowledgments 
and  this  explanation. 

Every  paper  is  preceded  by  an  abstract  in  a  standard  format.  Some  papers 
may  also  have  the  original  abstract  included.  Further,  you  may  recall  the 
Q&A  session  we  had  at  the  end  of  each  presentation.  In  cases  where  a 
technical  interchange  did  take  place,  the  questions  and  answers  are 
documented  at  the  end  of  each  pertinent  paper.  Every  paper  did  not  have 
a  Q&A  session,  and  I  have  included  all  Q&A  sheets  that  were  handed  to  me 
at  the  end  of  each  presentation.  Except  for  a  few  minor  editorial  changes, 
the  information  on  these  sheets  has  not  been  significantly  altered. 
Obviously,  these  sheets  are  as  "good"  as  the  inputs  you  provided. 

In  summary,  I  hope  the  above  explanation  was  helpful.  I  have  done  what  I 
consider  to  be  a  thorough  job  of  collecting  and  checking  all  the 
information  for  these  proceedings.  Errors  will  occur,  however,  and  while 
I  will  entertain  any  comments  and  criticisms  on  this  issue,  these 
proceedings  will  stand  as  published. 

Thank  you  for  your  participation,  and  your  patience! 


Capt  Vishnu  V.  Nevrekar  November  1987 

Earth  Sciences  Division 
USAF  Geophysics  Laboratory 
Hanscora  AFB,  MA  01731 


ACKNOWLEDGMENTS 


We  couldn't  possibly  organize  a  conference  the  scope  and  size  of  our  forum 
without  some  very  diligent  "behind-the-scenes"  work  by  a  few  outstanding 
individuals.  We  would  like  to  recognize  their  efforts  and  thank  them  for 
their  support  throughout  the  planning  and  execution  of  the  1987  Gradioraetry 
Conference. 

We  are  indebted  to  the  Directorate  of  Protocol  at  the  US  Air  Force  Academy 
for  allowing  us  to  host  the  Conference  there  these  past  13  years.  Ms  Nancy 
Gass  was  the  liaison  officer  from  the  Directorate,  and  we  gratefully 
acknowledge  her  assistance  in  handling  all  the  conference  arrangements, 
including  hotel  accomodation  for  the  attendees,  transportation  for  the 
Conference  and  NORAD/ CSOC  tours,  luncheons  during  the  conference  and  the 
"mixers"  later  in  the  evening,  all  of  which  were  set  up  with  great  skill 
and  professionalism. 

Also,  we  acknowledge  the  outstanding  efforts  of  TSgt  Kent  Droz,  USAF,  of  the 
Community  Relations  Division  at  HQ  NORAD,  who  arranged  tours  of  the  Cheyenne 
Mountain  Complex  (CMC)  and  the  Cbnsolidated  Space  Operations  Chnter  (CSOC) 
at  Falcon  AFS.  These  tours  gave  the  conference  attendees  a  first-hand  look 
at  the  complex  space  defense  environment. 

Next,  we  thank  all  the  speakers  for  taking  the  time  to  compile  and  present 
their  papers  for  the  benefit  of  the  (bnference  attendees.  As  always,  the  broad 
mix  of  topics  went  a  long  way  towards  making  the  Conference  an  intellectually 
stimulating  event.  Indeed,  the  high  quality  of  the  research  material  presented 
"made”  this  Conference. 

Finally,  we  thank  Col  J.R.  Johnson,  Commander,  AFGL,  Dr  Donald  H.  Eckhardt , 
Director,  Earth  Sciences  Division  and  Dr  Thomas  P.  Rooney,  Chief,  Geodesy 
and  Gravity  Branch,  without  whose  support  and  guidance  this  Conference  could 
not  have  been  held. 


Alphabetical  Listing  of  Conference  Participants 


Name 

*Georges  Balmino 
Anthony  Barringer 
*Hans  Baussus  Von  Luetzow 
Don  Benson 
Ed  Blegert 
John  Binns 
*Sam  Bose 
*Carl  Bowin 
John  Brozena 
Marcus  Chalona 
Lindrith  Cordell 
Ronald  Davis 
Mark  Dransfield 
Donald  Eckhardt 
Michael  Ellett 
Harry  Emrick 
John  Fett 
Charles  Finley 
Thomas  Flschetti 
James  Fix 
Cuy  Flanagan 
*Rene  Forsberg 
Caot  Terry  Fundak 
*David  Cleason 
Rob  Goldsborough 


Organization 

C.N.E.S. /Bureau  Gravimetrique  International  (FR) 

Barringer  Resources,  Inc 

US  Army  Engineer  Topographic  Laboratory 

Dynamics  Research  Corporation 

Shell  Development  Corporation 

BP  Minerals  International,  Ltd  (UK) 

Applied  Science  Analytics,  Inc 
Woods  Hole  Oceanographic  Institution 
Naval  Research  Laboratory 
US  Naval  Oceanographic  Office 
US  Geological  Survey 
Northrop  Electronics  Division 
University  of  Western  Australia  (AUS) 

USAF  Geophysics  Laboratory 
USAF  Space  Division 
Consultant 

LaCoste  and  Romberg  Gravity  Meters,  Inc 

National  Aeronautics  and  Space  Administration 

Technology  Management  Consultants,  Inc 

Teledyne  Geotech 

Standard  Oil  Production  Company 

Geodetic  Institute  (DEN) 

USAF  Geophysics  Laboratory 
USAF  Geophysics  Laboratory 
USAF  Geophysics  Laboratory 


Name 


Organization 


*John  Graham 
Andrew  Grierson 
Michael  Hadfield 
Richard  Hansen 
Chris  Harrison 
Ray  Hassanzadeh 
*Warren  Heller 
Howard  Heuberger 
George  Hinton 
Albert  Hsui 
Gene  Jackson 
Christopher  Jekeli 
♦Albert  Jircltano 
Col  J.R.  Johnson 
J.  Edward  Jones 
J.  Latimer 
Andrew  Lazarewlcz 
Thomas  Little 
♦Dan  Long 
James  Lowery  III 
Charles  Martin 
♦Bahram  Mashhoon 
♦Ernest  Metzger 
♦M.  Vol  Moody 
♦Ian  Moore 
ILt  Vishnu  Nevrekar 


Defense  Mapping  Agency 

Bell  Aerospace  Textron 

Honeywell,  Inc 

Colorado  School  of  Mines 

Geodynamics  Corporation 

McDonnell  Douglas 

The  Analytic  Sciences  Corporation 

Johns  Hopkins  University 

Consultant 

USAF  Geophysics  Laboratory 

McDonnell  Douglas 

USAF  Geophysics  Laboratory 

Bell  Aerospace  Textron 

Commander,  USAF  Geophysics  Laboratory 

USAF  Intelligence  Service 

Johns  Hopkins  University 

USAF  Geophysics  Laboratory 

US  Naval  Oceanographic  Office 

Eastern  Washington  University 

Rockwell  International 

University  of  Maryland  Research  Foundation 

University  of  Missouri 

Bell  Aerospace  Textron 

University  of  Maryland 

University  of  Queensland  (AUS) 

USAF  Geophysics  Laboratory 


♦Ho  Jung  Paik 


University  of  Maryland 


Name 

Maj  John  Prince 
*Rlchard  Rapp 
Richard  Reineman 
*Jean-Paul  Richard 
Thomas  Rooney 
Alan  Rufty 
Alton  Schultz 
*Michael  Sideris 
Ted  Sims 
Randall  Smith 
*David  Sonnabend 
Milton  Trageser 
Gary  Tuck 
Herbert  Valliant 
Robert  Valska 
*Frank  van  Kann 
Richard  Wold 
Robert  Ziegler 
*Alan  Zorn 
Paul  Zucker 


Organization 

USAF  Office  of  Scientific  Research 
Ohio  State  University 
GWR  Instruments 
University  of  Maryland 
USAF  Geophysics  Laboratory 
Naval  Surface  Weapons  Center 
AMOCO  Production  Company 
University  of  Calgary  (CAN) 

Naval  Surface  Weapons  Center 
Defense  Mapping  Agency 
CALTECH/ Jet  Propulsion  Laboratory 
Charles  Stark  Draper  Laboratory 
University  of  Queensland  (AUS) 

LaCoste  and  Romberg  Gravity  Meters,  Inc 

Defense  Mapping  Agency 

University  of  Western  Australia  (AUS) 

TerraSense,  Inc 

Defense  Mapping  Agency 

Dynamics  Research  Corporation 

Johns  Hopkins  University 


*  Denotes  Speaker  at  Conference 


OVERALL  DATA  PROCESSING  SCHEME 


APPLIED  SCIENCE  ANALYTICS.  INC. 


GRADIOMETER  INLINE  AND  CROSSLINE  MEASUREMENTS 


APPLIED  SCIENCE  ANALYTICS,  INC. 


GRAVITY  GRADIOMETER  DATA  PROCESSING 
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Calibration  Overview 

LABORATORY  SELF  GRADIENT  DATA  CONSISTING  OF  AVERAGED  6GI  OUTPUTS  TAKEN 
OVER  280  COMBINATIONS  OF  SCORSBY  AZIMUTH.  PITCH  AND  CAROUSEL  ANGLES. 

LABORATORY  STEP  INDEXING  AND  CONTINUOUS  CAROUSELLEO  DATA. 

VEHICLE  SELF  GRADIENT  DATA  CONSISTING  OF  AVERAGED  GGI  OUTPUTS  TAKEN  OVER 
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GGI  #1  OUTPUT  DURING  U\BORATORY  CALIBRATION  WITH  OUTPUT 
PREDICTED  BY  MASS  HOOEL  TAYLOR  SERIES 


612.5 


INLINE  I 


9.' J  J 


CROSS  1 


GGI  #2  OUTPUT  DURING  LABORATORY  CALIBRATION  WITH  OUTPUT 
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GGI  #1  OUTPUT  DURING  AIRCRAFT  CALIBRATION  WITH  OUTPUT 
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TECHNICAL  APPROACH 
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TITLE  OF  PAPER;  Gravity  Gradiometer  Survey  System  (GGSS)  Data  Processing 
and  Data  Use 

SPEAKER;  Warren  G.  Heller 
QUESTIONS  AND  COMMENTS; 


1 .  Question;  Unknown 

How  big  is  the  test  area  you  used? 

Response; 

300  km  on  a  side,  modelled  over  40  km. 

2.  Question;  Rene  Forsberg 

In  fitting  your  improved  AWN  covariance  model,  did  you  use  local  gravity 
data  (wavelengths  shorter  than  10  km)  to  compare  "topographic"  and 
"non-topographic"  local  gravity  power? 

Response ; 

No,  these  data  were  not  available.  However,  the  slope  of  the  refined  model 
seems  to  fit  well  with  the  slope  inferred  by  the  original  AWN  model. 

Question;  Richard  Rapp 

What  is  the  accuracy  of  the  recovery  of  the  terrain  signal  in  the  block 
sizes  you  considered  to  be  the  resolution  of  the  system? 

Response ; 

The  aim  is  to  recover  the  terrain  effects  to  0.1  mgal. 

4.  Question;  Chris  Jekeli 

Did  you  use  isostatic  compensation  model  for  computing  terrain  effects  on 
deflection  of  vertical? 

Response; 

No,  just  used  the  terrain  data. 

5.  Question;  John  Brozena 

Did  the  error  model  for  the  gradiometer  used  in  your  analysis  include 
environmental  noise  sources? 

Response ; 


Yes  • 


6 .  Question : 


J  ame  s  E •  Fix 


In  analyzing  the  terrain  effect,  was  a  variable  density  or  a  constant  density 
used? 

Response : 

A  constant  density  was  used.  Density  was  taken  as  2.67  g/cm^. 

7.  Question :  Sam  Bose 

How  does  the  performance  deteriorate  as  the  averaging  size  is  increased? 
Response : 

There  is  no  performance  degradation  provided  the  averaging  size  takes  into 
account  the  maximum  bandlimit  of  the  gradiometer  signal. 


SELF  GRADIENT  CALIBRATION  OF  THE  GGSS 
IN  A  C-130  AIRCRAFT 


by 

Dr.  W.  John  Hutcheson 
Bell  Aerospace  Textron 
P.O.  Box  One 
Buffalo,  NT  142A0 


ABSTRACT 

Due  to  the  inverse  cube  law  for  gravity  gradients,  mass  structures  close 
to  the  gradlometer  sensing  elements  produce  significant  outputs  termed  self 
gradients  which  have  to  be  compensated  for  in  the  GGSS  Stage  I  data  reduction. 

In  the  Bell  approach  to  the  self  gradient  calibration,  a  mass  model  representing 
the  mass  structures,  consisting  of  the  GGSS  platform  gimbals,  servo  motors, 
binnacle  and  aircraft,  la  identified  using  optimal  identification  techniques 
and  then  used  to  generate  the  compensation. 

This  paper  contains  a  brief  description  of  the  theory  underlying  the  Bell 
approach  to  the  self  gradient  calibration,  details  of  the  self  gradient  calibra¬ 
tion  of  the  GGSS,  covariance  results,  GGSS  calibration  data  and  the  calibration 
curves  representing  the  combined  field  of  the  GGSS,  the  van  and  the  C-130 


aircraft . 


MOVING  BAS 


to  temperature  fluctuations.  Improvements  to  the  cryostat  have  increased  the  thermal  isolation 
and  stability  and  are  expected  to  reduce  noise  from  this  source. 


Figure  5  Gradient  signal  as  a  function  of  distance.  The  observed  gradient  signal  is  plotted 
against  the  distance  between  the  centres  of  mass  of  the  gradiometer  and  the  rim  of  the 
gradient  generator.  The  solid  curve  is  calculated  from  its  known  mass  distribution 
and  varies  approximately  as  the  inverse  cube  of  the  spacing. 


Conclusion  ' 

The  prototype  gradiometer  has  demonstrated  that  useful  gradient  sensitivity  is  attainable  and 
that  intrinsic  detector  noise  is  unlikely  to  be  a  limitation  on  the  development  of  a  practical 
instrument.  However,  improvements  in  multi-axis  common  mode  rejection  and  rotational 
stabilization  are  required.  These  will  be  difficult  to  achieve,  but  do  appear  to  be  practicable. 
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TITLE  OF  PAPER:  A  Prototype  Superconducting  Gravity  Gradiometer 
for  Geophysical  Exploration 


SPEAKER:  Frank  J.  van  Kann 

QUESTIONS  AND  COMMENTS: 

1.  Question:  Jean-Paul  Richard 

What  are  the  mechanical  Q,  the  resonant  frequency  of  test 
masses  and  the  electrical  Q? 

Response : 

^mechanical  " 

Frequency  of  test  masses  =  30  Hz 
Qelectrical  ~  not  determined 
2-  Question:  Ho  Jung  Paik 


a)  What  kind  of  suspension  did  you  use  for  the  gradiometer? 

b)  Why  Is  It  that  your  noise  spectrum  does  not  show  resonance  peaks  of  the 
suspension  nodes? 


Response: 


a)  It  was  soft-suspended.  And  there  was  some  nasty  resonance. 

b)  We  used  notch  filters  to  remove  those  peaks. 


3.  Question;  Anthony  R.  Barringer 


What  type  of  rotational  stability  do  you  require  in  the  platforn  for 
mounting  the  gradiometer? 


Response : 


10“^  radians  per  second.  It  is  desirable  to  have  the  platform  inside  the 
cryostat . 


GRAVITY  GRADIOMETER  SURVEY  SYSTEM  (GGSS) 
DATA  PROCESSING  AND  DATA  USE 


by 

Dr.  Warren  G.  Heller 
The  Analytic  Sciences  Corp. 
100  Walkers  Brook  Drive 
Reading,  MA  01867 


ABSTRACT 

Since  the  GGSS  will  be  flown  at  a  given  altitude,  h,  (approx.  600tn)  above 
the  surface,  a  short  wavelength  limit  is  effectively  imposed  on  the  information 
content  of  the  acquired  data.  This  limit  Is  dictated  by  the  noise  of  the 
gradiometer  instruments  and  the  upward  continuation  factor,  where  X  is 

the  gravity  disturbance  wavelength.  Since  the  information  is  band  limited,  it 
is  appropriate  to  consider  representing  the  downward  continued  gravity  disturb¬ 
ance  estimates  as  area  means  over  a  suitably-sized  block  that  retains  full  data 
resolution  and  is  easy  to  incorporate  into  existing  gravity  data  bases.  For  a 

given  survey  area,  the  averaging  block  size  increases  with  flight  altitude. 

This  paper  1)  describes  an  analytic  technique  for  determining  the  shortest  wave¬ 
length  at  which  information  is  reliably  gathered  by  an  airborne  gradiometer,  2) 
presents  the  results  of  applying  this  technique  in  the  GGSS  test  area,  and  3) 
discusses  the  implications  of  survey  altitude  on  resolution  of  gravity  disturb¬ 
ance  recovery  by  gradiometrlc  surveys  in  other  areas.  Video  displays  are 
presented  which  Illustrate  character  of  the  short  wavelength  gravity  field  in 


the  test  area. 
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FOREWORD 


This  document  contains  material  used  in  a 
presentation  given  by  The  Analytic  Sciences  Corpo¬ 
ration.  The  material  is  not  intended  to  be  self- 
explanatory,  but  rather  should  be  considered  in  the 
context  of  the  overall  presentation. 


THE  ANALYTIC  SCIENCES  CORPOnATION 


ABSTRACT 


Gravity  Gradiometer  Survey  System  (GGSS) 
Data  Processing  and  Data  Use 


Since  the  GGSS  will  be  flown  at  a  given  altitude,  h, 
(approx.  600  m)  above  the  surface,  a  short  wavelength  limit  is 
effectively  imposed  on  the  information  content  of  the  acquired 
data.  This  limit  is  dictated  by  the  noise  of  the  gradiometer 

instruments  and  the  upward  continuation  factor,  e  where 
\  is  gravity  disturbance  wavelength.  Since  the  information  is 
band  limited,  it  is  appropriate  to  consider  representing  the 
downward  continued  gravity  disturbance  estimates  as  area  means 
over  a  suitably-sized  block  that  retains  full  data  resolution  and 
is  easy  to  incorporate  into  existing  gravity  data  bases.  For  a 
given  survey  area,  the  averaging  block  size  increases  with  flight 
altitude.  This  paper  1)  describes  an  analytic  technique  for 
determining  the  shortest  wavelength  at  which  information  is  reliabl 
gathered  by  an  airborne  gradiometer,  2)  presents  the  results 
of  applying  this  technique  in  the  GGSS  test  area,  and  3)  discusses 
the  implications  of  survey  altitude  on  resolution  of  gravity 
disturbance  recovery  by  gradiometric  surveys  in  other  areas. 

Video  displays  are  presented  which  illustrate  the  character  of 
the  short  wavelength  gravity  field  in  the  test  area. 
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THERMAL  NOISE  IS  ONLY  A  SMALL  CONTRIBUTOR  TO 
TERRESTIAL  MOVING  BASE  GRAVITY  GRADIOMETERS 


TITLE  OF  PAPER;  Bell  Aerospace  Gravity  Gradlometer  Survey  System  (GGSS)  - 
Program  Review 

SPEAKER:  Ernest  H.  Metzger 

QUESTIONS  AND  C0^C1E^TS: 

1.  Question:  David  Gleason 

Does  the  van  have  a  cruise  control  system? 

Response : 

Yes,  but  the  cruise  control  system  will  be  of  assistance  only  if  the  van's 
velocity  >  20  mph. 

2.  Question:  Anthony  R.  Barringer 

How  many  satellites  did  you  use  for  GPS  positioning? 

Was  a  ground  reference  station  used? 

Response: 


Four  satellites  as  a  minimum  essential.  No  ground  receiver  used. 

3.  Question:  Jim  Lowery 


What  was  used  as  an  altitude  reference  on  the  results  shown  crossing  the 
gorge? 


Response: 


GPS  aided  by  barometric  altimeter  was  used  as  the  altitude 
reference  with  an  accuracy  of  7  meters  rms  to  20  meters  rms. 


4.  Question:  Ted  Sims 


What  fraction  of  a  ”g"  can  you  reasonably  expect  the  system  to 
experience  while  in  a  turn? 


Response :  .5  g  is  reasonable. 


A  PROTOTYPE  SUPERCONDUCTING  GRAVITY  GRADIOMETER 
FOR  GEOPHYSICAL  EXPLORATION 


by 

Dr.  Frank  J.  van  Kann,  et  al 
University  of  Western  Australia 
Department  of  Physics 
Nedlands,  Western  Australia  6009 


ABSTRACT 

A  three  axis  gradiometer,  designed  to  measure  the  diagonal  components  of 
the  earth's  gravitational  gradient  tensor,  has  been  built  and  is  being  tested 
in  the  laboratory.  It  consists  of  three  pairs  of  accelerometers.  The  acceler¬ 
ometers  of  each  pair  are  mounted  with  their  sensitive  axes  co-linear  and 
orthogonal  to  the  other  pairs.  The  difference  in  acceleration  for  a  pair  is 
proportional  to  the  appropriate  component  of  the  gradient  tensor  and  is  sensed 
via  a  displacement  which  modulates  the  inductance  of  a  superconducting  coil 
coupled  by  means  of  a  transformer  to  an  RF  biased  SQUID  with  energy  sensitivity 
3  X  10-29  J/Hz. 

Rejection  of  in-line  common  mode  acceleration  is  achieved  by  tuning  the 
natural  resonant  frequencies  of  the  accelerometers  by  adjustment  of  persistent 

currents  stored  in  the  superconducting  force  coils.  A  common  mode  rejection  ratio 
near  100  dB  has  been  achieved  in  the  presence  of  common  mode  accelerations 
approaching  10”^  ms“^.  This  has  enabled  the  detection  of  a  laboratory  generated 
signal  as  small  as  5  Ed  at  signal  frequencies  below  1  Hz  with  signal  to 
noise  ratio  approaching  10.  Above  0.1  Hz,  the  noise  floor  of  the  instrument  is 
about  0.5  Ed  //  Hz  under  quiet  conditions.  Below  0.1  Hz  it  has  been  limited 

by  thermal  drifts  but  measurements  are  at  present  being  carried  out  in  a  new 
cryostat  with  improved  temperature  stability. 


A  PROTOTYPE  SUPERCONDUCTING  GRAVITY  GRADIOMETER 


FOR  GEOPHYSICAL  EXPLORATION* 


F  J  van  Kann,  M  J  Buckingham,  M  H  Dransfield,  C  Edwards,  A  G  Mann,  R  D  Penny  and  P  J  Turner 
Physics  Department,  The  University  of  Western  Australia,  Nedlands,  6009,  Australia. 


Abstract 

A  three  axis  gradiometer,  designed  to  measure  the  diagonal  components  of  the  earth's 
gravitational  gradient  tensor,  has  been  built  and  is  being  tested  in  the  laboratory.  It  consists  of 
three  pairs  of  accelerometers.  The  accelerometers  of  each  pair  are  mounted  with  their  sensitive 
axes  co-linear  and  orthogonal  to  the  other  pairs.  The  difference  in  acceleration  for  a  pair  is 
proportional  to  the  appropriate  component  of  the  gradient  tensor  and  is  sensed  via  a 
displacement  which  modulations  the  inductance  of  a  superconductino  coil  coupled  by  means  of  a 
transformer  to  an  RF  biased  SQUID  with  energy  sensitivity  3  x  10'*®  J/Hz. 

Rejection  of  in-line  common  mode  acceleration  is  achieved  by  tuning  the  natural  resonant 
frequencies  of  the  accelerometers  by  adjustment  of  persistent  currents  stored  in 
superconducting  force  coils.  A  common  mode  rejection  ratio  near  1 00  dB  has  been  achieved  in 
the  presence  of  common  mode  accelerations  approaching  10*^  ms’^.  This  has  enabled  the 
detection  of  a  laboratory  generated  signal  as  small  as  5  Ed  at  signal  frequencies  below  i  Hz  with 
signal  to  noise  ratio  approaching  10.  Above  0.1  Hz,  the  noise  floor  of  the  instrument  is  about 
0.5  Ed/VHz  under  quiet  conditions.  Below  0.1  Hz  it  has  been  limited  by  thermal  drifts  but 
measurements  are  at  present  being  carried  out  in  a  new  cryostat  with  improved  temperature 
stability. 

Introduction 

The  form  of  the  earth's  gravitational  potential  function  contains  a  wealth  of  information  of 
importance  in  geophysics.  For  the  purposes  of  geophysical  exploration,  this  has  traditionally 
been  exploited  through  measurement  of  the  first  spatial  derivatives  of  the  potential  -  the 
gravity  field.  Because  of  the  difficulty  of  distinguishing  spatial  variations  in  gravity  from 
temporal  fluctuations  of  the  acceleration  of  a  moving  vehicle,  these  measurements  of  gravity  can 
be  made  to  sufficient  precision  only  with  stationary,  earth  based  instruments.  The  limitations 
imposed  by  translational  acceleration  can  in  principle  be  avoided  by  measurement  of  the  second 
derivative  of  the  potential  •  gravity  gradients.  Indeed,  the  discrimination  of  interestina 
geological  anomalies  could  be  more  easily  achieved  by  direct  measurements  of  the  gradient^ 
rather  than  gravity  itself  and  under  appropriate  conditions  gradient  measurements  are  less 
dependent  on  elaborate  corrections  for  topographical  features^. 

To  obtain  gravity  gradient  data  useful  for  exploration,  a  noise  level  less  than  0.1  feo/\'Hz  is 
required,  which  implies  an  equivalent  acceleration  resolution  on  the  order  of  10'^^  ms'^  in  an 
instrument  of  reasonable  size  and  mass.  The  extremely  large  common  mode  rejection  ratio 
(possibly  exceeding  200  dB)  required  to  make  these  measurements  in  a  moving  vehicle  may  be 
attainable,  given  a  system  with  adequate  linearity  and  a  sufficiently  precise  and  stable  method  of 
tuning.  However,  the  finite  elastic  stiffness  of  materials  gives  rise  to  errors  in  the  gradient 


signal  which  are  quadratic  in  the  common  mode  acceleration.  The  size  of  these  errors  depends  on 
the  geometrical  shape  and  the  elastic  stiffness  of  the  instrument,  but  for  materials  with  a 
velocity  of  sound  around  3  km/s  and  reasonable  shape  the  maximum  allowable  common  mode 
acceleration  is  less  than  10'^  ms*^  for  a  0.1  Ed  error.  This  sets  an  exacting  requirement  for 
the  translational  acceleration  isolation  of  the  stabilisation  system  required  for  any  vehicle 
suitable  for  use  as  an  exploration  platform. 

Rotational  stabilisation  is  also  required  to  reduce  errors  which,  for  the  diagonal  components  of 
the  gradient  tensor,  are  quadratic  in  the  angular  velocity  of  the  instrument.  For  these  errors  to 
be  less  than  0.1  E6  requires  the  angular  rate  to  be  less  than  10'®  radian  s'^  about  any  axis. 
Rotation  sensors  with  adequate  performance  to  meet  this  requirement  are  currently  available. 
However,  these  will  need  to  be  adapted  for  low  temperature  operation,  since  the  innate 
mechanical  elastic  compliance  of  the  cryostat  imposed  by  thermal  design  considerations  will 
require  that  some  rotational  stabilisation  be  Implemented  inside  the  cryogenic  environment. 

The  laboratory  prototype 

The  three  axis  prototype  gradiometer  uses  principles  similar  to  those  described  by  Mapoles^, 
Paik'^  and  Moody  et  aP.  It  consists  of  six  essentially  identical  accelerometers  grouped  to  form 
three  pairs,  one  for  each  tensor  component  to  be  sens^.  The  two  end  faces  of  each  accelerometer 
can  be  identified  by  a  letter  A,  B,  C  or  0  so  that  the  two  accelerometers  for  a  given  pair  can  be 
labelled  AB  and  CD  respectively.  These  are  selected  for  matched  mechanical  resonant  frequencies 
and  are  mounted  with  their  sensitive  axes  co-llnear  and  orthogonal  to  those  of  the  remaining 
pairs.  Each  accelerometer  consists  of  a  solid  niobium  cylinder,  some  30  mm  in  diameter,  30 
mm  in  length,  and  about  300  gm  in  weight  suspended  at  each  end  by  a  thin  folded  cantilever 
niobium  leaf  spring  in  a  niobium  housing. 

The  remote  end  faces  A  and  0  of  the  pair  are  parallel  and  in  close  proximity  to  annular,  single 
layer,  spiral  "pancake"  niobium  wire  coils  attached  to  the  ends  of  the  housing.  Each  of  these 
pancake  coils  actually  consists  of  a  pair  of  concentric  coils:  the  smaller  inner  one  being  used  for 
RF  position  sensing  and  the  larger  outer  one  forming  the  force  coils  for  CMRR  tuning  and 
feedback.  The  resonant  frequency  of  the  accelerometers  is  about  25  Hz  and  can  be  increased  by 
several  percent  by  means  of  a  persistent  current  stored  in  the  appropriate  force  coil.  The  end 
face  labelled  C  of  accelerometer  CO  is  similarly  with  another  pancake  coil  mounted  on  the  B  end 
face  of  the  paired  accelerometer  AB.  This  coil  is  coupled  to  the  SQUID  by  means  of  a 
superconducting  matching  transformer  and  Is  used  to  directly  sense  the  differential  motion 
between  the  accelerometers. 

Common  mode  acceleration  sensing 

Accelerations  are  monitored  by  sense  coils  at  the  ends  of  the  gradiometer  housing.  Each  sense 
coil  is  incorporated  into  the  tank  circuit  of  a  radio  frequency  oscillator,  whose  frequency  is 
modulated  by  motion  of  the  test  mass  relative  to  the  housing.  This  position  readout  permits 
preliminary  testing  of  the  accelerometers  at  room  temperature  and  also  enables  calibration  of 
the  primary  superconducting  differential  motion  sensing  system  when  cold. 

The  RF  oscillators  have  been  optimised  for  low  power  operation,  both  to  permit  their  use  in  the 
highly  thermally  isolated  cryogenic  environment  and  also  to  minimise  SQUID  interference.  At 
liquid  helium  temperature,  these  provide  a  stable  and  sensitive  position  sensor  with  10'^  m 


resolution  while  dissipating  only  40  pW. 


The  spring  constant  of  the  mechanical  springs  is  augmented  as  required  by  means  of  the  magnetic 
force  from  a  persistent  current  stored  in  the  force  coils.  This  allows  the  accelerometers  to  be 
precisely  matched  to  achieve  high  rejection  of  accelerations  along  the  gradient  sensing  axis. 

The  effectiveness  of  the  adjustment  of  persistent  current  for  CMRR  tuning  is  illustrated  in 
figure  1 .  Here  the  natural  resonant  frequency'of  the  CD  accelerometer  was  about  0.8  Hz  below 
that  of  AB.  The  upper  two  curves  in  figure  1  show  the  Fourier  spectrum  of  the  response  of  the 
two  accelerometers  from  white  noise  excitation.  The  peak  near  26.5  Hz  corrresponds  to  the 
natural  untuned  low  frequency  normal  mode  for  the  coupled  oscillators.  (The  other,  high 
frequency  normal  mode  is  above  28  Hz  and  not  visible  in  this  diagram.)  The  family  of  curves  in 
the  central  region  of  figure  1  shows  the  spectrum  of  the  transfer  function  amplitude  i.e.  the 
magnitude  of  the  complex  ratio  of  the  response  of  the  two  accelerometers,  for  several  persistent 
currents  stored  in  force  coils  C  and  D.  The  lower  curves  show  the  corresponding  phase  of  the 
transfer  function. 


Figure  1  Fourier  spectrum  of  the  gradiometer  response  as  a  function  of  tuning  current. 


The  curves  1  to  5  show  the  resonant  frequency  of  the  CD  accelerometer,  corresponding  to  tho 
peak  in  the  transfer  function  amplitude,  being  increased  to  match  that  of  AS.  corresponding  to 
the  dip  in  amplitude.  The  accelerometers  are  matched  when  these  coincide  as  in  curve  number  4. 
In  curves  1  to  3,  the  stored  current  is  too  small  and  the  frequency  of  CD  is  lower  than  that  of  AB. 
In  curve  5.  the  current  is  loo  large  and  the  frequency  difference  is  reversed.  In  curve  4,  the 
accelerometers  are  as  closely  matched  as  can  be  determined  by  this  method.  More  precise  tuning 
is  achieved  by  direct  measurement  of  the  differential  motion  using  the  SQUID. 

Differential  mode  acceleration  sensing 

Differential  motion  between  the  two  accelerometers  is  measured  to  extremely  high  resolution  by 
means  of  an  RF  biased  SQUID  magnetometer,  model  330X.  manufactured  by  Biomagnetic 
Technology  Inc.  This  detects  changes  in  the  persistent  current  trapped  in  the  superconducting 
differential  motion  sense  coil.  Extreme  care  has  been  taken  in  shielding  this  input  circuitry 
from  fluctuations  in  the  ambient  magnetic  field  and  also  from  RF  interference  which  can  cause 
the  SQUID  to  cease  functioning. 


Figure  2  Gradiometer  response  to  applied  calibration  and  common  mode  rejection  test  signals. 


The  differential  motion  sensitivity  is  calibrated  with  respect  to  the  known  sensitivity  of  the  RF 


position  sensor  as  shown  in  figure  2.  The  Fourier  components  of  both  the  common  mode  motion, 
sensed  by  the  RF  detector  and  differential  mode  motion,  sensed  by  the  SQUID  are  shown  as  a 
function  of  frequency.  The  vertical  scale  on  the  left  is  labelled  with  the  calibrated  displacement 
scale,  while  the  scale  on  the  right  shows  the  equivalent  acceleration  amplitude  relative  to  the 
earth's  gravitational  acceleration.  For  the  calibration,  the  CD  accelerometer  is  forced  into 
oscillation  at  a  known  amplitude  and  frequency  by  means  of  one  of  its  force  coils,  which  has  a 
trapped  persistent  current.  Although  the  heat  switch  which  is  in  parallel  with  the  force  coil 
remains  cold,  the  stray  series  inductance  allows  the  persistent  current  to  be  modulated  via  Iho 
external  current  leads.  Since  the  forced  oscillation  is  well  below  the  resonant  frequency,  the 
resulting  motion  of  the  AB  accelerometer  is  small  and  not  detectable  above  the  noise.  The  CMRR 
is  measured  simultaneously  by  means  of  a  forced  common  mode  oscillation  of  the  entire 
gradiometer  assembly,  which  itself  is  suspended  inside  the  dewar  on  soft  coil  springs  with  a 
resonant  frequency  of  about  1  Hz.  The  dependence  of  the  CMRR  on  the  trapped  current  is  shown 
in  figures.  The  circles  and  crosses  represent  data  from  two  different  runs,  with  Iq^  =  loo 

and  20  A^  respectively.  For  clarity,  some  of  the  data  from  the  latter  are  omitted  and  plotted  on 
an  expanded  scale  in  the  inset.  The  maximum  CMRR  achieved  is  nearly  100  dB. 


-100  0  100  200 


oeVIATION  FROM  OPTIMUM 
TUNING  CURRENT. 

Figure  3  Common  mode  rejection  ratio  as  a  function  of  the  square  of  the  push  coil  current. 


A  gravitational  gradient  generator  was  constructed  to  test  the  performance  of  the  gradiometer 
with  actual  time  varying  gradients.  The  generator  consists  of  a  1 .2  meter  diameter  wheel  at  the 
periphery  of  which  are  attached  four  lead  masses  weighing  some  65  kg  each.  The  gradient 
produced  by  such  a  mass  when  placed  close  (0.3  m)  to  the  gradiometer  is  some  120  EO.  When 
the  disc  is  set  into  rotation  by  a  variable  speed  drive,  it  produces  an  AC  gradient  with 


fundamental  Fourier  component  at  four  times  the  rotation  frequency,  and  rms  amplitude  of  about 
30  Ed. 


The  ability  of  the  gradiometer  to  successfully  detect  the  gravitational  gradient  of  the  generator  is 
shown  clearly  in  figure  4.  The  frequency  for  the  measurement  was  chosen  such  that  the 
fundamental  rotor  frequency  (0.077  Hz)  and  its  first  few  harmonics  did  not  coincide  with  any 
natural  resonances  of  the  gradiometer  suspension  or  dewar  system.  The  gradient  signal  at 
0.3  Hz  has  the  expected  amplitude  of  approximately  4  E6  rms.  This  fundamental  Fourier 
component  of  the  gradient  produced  by  the  rotor  can  be  quantitatively  predicted  from  its  known 
mass  distribution,  and  serves  as  a  useful  check  on  the  gradiometer  calibration.  The  strong 
signal  at  0.53  Hz  results  from  the  rocking  motion  of  the  gradiometer  on  its  suspension  springs 
and  cannot  be  suppressed  by  common  mode  rejection. 


Figure  4  Detection  of  a  gradient  signal.  This  figure  shows  on  a  logarithmic  scale  the  Fourier 
spectrum  of  the  differential  motion  SQUID  signal  in  response  to  the  four  mass 
gradient  generator  rotating  at  0.077  Hz.  The  SQUID  sensitivity  is  2  V/^iA  and  the 
unfiltered  output  is  used  for  feedback  damping  to  reduce  the  accelerometer  Q.  The 
signal  is  processed  by  an  HP  3582A  spectrum  analyser  without  any  additional 
filtering. 


The  distance  dependence  of  the  gravity  gradient  produced  by  the  generator  (approximately  i/r  ' 
at  short  distances)  is  easily  calculated  and  the  results  are  in  good  agreement  with  the  observed 
behaviour,  as  shown  in  figure  5. 

At  very  low  frequency,  the  gradiometer  noise  level  rises  significantly  because  of  the  sensitivity 
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Normal  flotation  level  5 


Flotation  level  with  spring  5’ 


Required  spring  constant  <0.14  Nrh 


TITLE  OF  PAPER:  A  Mercury  Manorneter  Gravity  Gradiometer 
SPEAKER:  G.  Ian  Moore 

OIT.STIOMS  AND  CO^^^fE^:TS: 

1.  Quest  ion :  Jean-Paul  Richard 

Sensitivity  of  capacitance  detection? 

Response : 

10~^  times  the  gap  =  2  x  10“^ 


mm. 


BELL  AEROSPACE  GRAVITY  GRADIOMETER  SURVEY  SYSTEM  (GGSS) 

PROGRAM  REVIEV’J 


by 


Mr.  Ernest  H.  Metzger 
Mr.  Louis  L.  Pfohl 

Bell  Aerospace  Textron 
P.O.  Box  One 
Buffalo,  MY  14240 


ABSTRACT 

A  review  of  GGSS  program  activities  in  1986  includes  system  lab  testing, 
land  vehicle  and  aircraft  installations,  electrical  power  and  signal  interfacing, 
and  shakedown  cruises.  Among  the  significant  accomplishments  were  system  output 
noise  determination  in  the  laboratory,  platform  and  aircraft  self-gradient 
calibrations,  and  Implementation  of  automated  flight  pattern  control  via  GGSS 
navigator  and  computer  linked  to  the  C-130  autopilot. 


MOVING  BASE 
GRAVITY  GRADIOMETER 


Bell  Aerospace 


986  GGSS  Activities 


Bell  Aerospace 


Functional  Block  Diagram 


SOLUTION  OF  MEASUREMENT  INTEGRALS  WITHOUT  MATRIX  INVERSIONS 
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GRAVITY  VECTOR  EXPRESSED  IN  TERMS  OF  K-L  COEFFICIENT  ESTIMATES 


APPLIED  SCIENCE  ANALYTICS,  INC. 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  Of  STANDARDS- 1963-A 


INTERPOLATED  AND  DOWNWARD  CONTINUED  2D  GRID  ESTIMATES  OF  GRAVITY  VECTOR 
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TttLE  of  PAPER:  Gravity  Gradioraeter  Data  Processing  Using 
the  Karhunen-Loeve  Method 

SPEAKER :  Sa;n  C.  Bose 

nUESTIOVS  A>n)  COM>fEXTS: 

1.  Question :  Ted  Sims 

Will  the  method  presented  accomodate  data  taken  at  differing  altitudes? 
Response : 

Yes. 

2 .  Question :  Hans  Baussus  von  Luetzow 

How  do  you  cousider  gravity  gradiometer  red  noise? 

Response : 

Red  noise  and  white  noise  effects  are  integrated  in  one  error  variance. 

3.  Question .  Anthony  R.  Barringer 

I  am  not  clear  on  your  survey  pattern  of  flying.  Do  you  have  a  viewgraph 
Response : 


Mo.  My  analysis  is  based  on  an  orthogonal  grid. 


NUMERICALLY  DERIVING  THE  KERNELS  OF  AN  INTEGRAL  PREDICTOR 
YIELDING  SURFACE  GRAVITY  DISTURBANCE  COMPONENTS 
FROM  AIRBORNE  GRADIENT  DATA 


by 

David  M.  Gleason 
Geodesy  and  Gravity  Branch 
Earth  Sciences  Division 
Air  Fbrce  Geophysics  Laboratory 
Hanscom  AFB,  MA  01731-5000 


ABSTRACT 

C.  Jekeli  (1986)  developed  an  integral  estimator  which,  when  used  in 
conjunction  with  a  set  of  airborne  gradient  observations,  yields  gravity  distur¬ 
bance  component  differences  between  a  desired  collection  of  actual  disturbance 
component  values  on  the  ground  and  a  corresponding  collection  of  least-squares 
collocation  predicted  values  that  are  based  on  a  small,  given  set  of  disturbance 
component  tie  point  values,  also  on  the  ground,  which  provide  needed  long-wave¬ 
length  gravity  Information.  (Thus  the  desired  actual  values  can  be  estimated 
by  adding  back  the  differences  to  the  tie  point-implied  values).  This  paper 
shows  how  all  18  possible  kernels  of  the  integral  estimator  can  be  easily  and 
accurately  approximated  via  two  dimensional  discrete  inverse  Fourier  transforms. 
Armed  with  such  a  set  of  kernel  values,  a  few  tie  points  and  a  set  of  airborne 
gradient  values  implied  by  a  mass  layer  gravity  model  for  northern  Texas,  the 
RMS  error  of  a  set  of  predicted  ground  disturbance  components,  referenced  to 
"'’rue"  values  Implied  by  the  same  model,  is  less  than  1  mgal.  A  flat  earth 
approximation  is  employed  in  this  exercise  using  (X,  Y,  Z)  (east,  north,  and 


down)  coordinates. 


NUMERICALLY  DERIVING  THE  KERNELS  OF  AN 
INTEGRAL  PREDICTOR  YIELDING  SURFACE  GRAVITY 
DISTURBANCE  COMPONENTS  FROM  AIRBORNE 


GRADIENT  DATA. 


•  FOR  A  DETAILED  EXPLANATION  OF  THE  NSWC  MASS  LAYER  LOCAL 
GRAVITY  MODEL  FOR  NORTHERN  TEXAS,  SEE  WHITE  (1984)  • 

(TASC/AFGL-TR-85-0037) _ 

I.  GETTING  INITIAL  LSC/TIE  POINT  PREDICTIONS: 

IF  U  CONTAINS  A  FEW  GIVEN  TIE  POINT  DISTURBANCE  COMPONENT 
VALUES  ON  THE  GROUND,  I.E., 

U  =  Ij  (J  =  X,YORZ) 

(Nj  by  I ) 

WE  CAN  ALWAYS  PREDICT  A  VECTOR  W  OF  OTHER  GROUND 
DISTURBANCE  COMPONENTS  VIA  THE  LSC  EQUATION 

W  =  [PI U  (  I  ) 

(Np  BY  I ) 

WHERE  THE  Np  BY  Nj  ESTIMATOR  MATRIX 

[PI  =  [C^l  ■  [  [Cyyl^-[0]1-'  <2) 

it 


(  Np  BY  Nj  ) 


(Nj  by  Nj) 


II  ESTIMATING  THE  (ACTUAL-LSC^  DISTURBANCE  COMPONENT 
DIFFERENCES 

C.  JEKELI  ( 1 986)  SHOWS  IE 

(1)  THE  VECTOR  y  CONTAINS  A  SET  OF  FAIRLY  DENSE  AND 
INFINITELY  EXTENDED  AIRBORNE  GRADIENTS  AT  SOME  CONSTANT 
ALTITUDE  H  ABOVE  THE  PLANE  EARTH  (THE  OBSERVATIONS  MAY 
OR  MAY  NOT  BE  REGULARLY  GRIDDED) 

(2)  THE  VECTOR  U  CONTAINS  A  FEW  TIE  POINT  DISTURBANCE 
COMPONENT  VALUES  ON  THE  GROUND  AT  POINTS  (Xk,Y|(,0),K=  I  ,Nj 

(3)  WE  ASSUME  THE  GRAVITY  SIGNAL  TO  BE  STATIONARY  AND  THE 
COVARIANCES  TO  BE  FUNCTIONS  ONLY  OF  THE  DISTANCE  BETWEEN 
A  PAIR  OF  POINTS  (P|  ,P2) 


THEN  THE  DIFFERENCES  BETWEEN  THE  ACTUAL  DISTURBANCE 


COMPONENT  VALUES  ON  THE  GROUND  AND  THE  CORRESPONDING 
VALUES  IMPLIED  BY  THE  LSC  PREDICTOR  OF  SECTION  I ,  AT  THE 
PREDICTION  POINT  (Xq,Yq,  0),  CAN  BE  EXPRESSED  AS 


00  oo 


W(XQ,Yp,0)  =  I  j[B(XQ-X,YQ-V,0)] •V(X,Y,H)dXdY  ^ 


•-O0  •■OO 


(  3  ) 


WHERE  EACH  INDIVIDUAL  KERNEL  ELEhENT  IN  THE  MATRIX  [  B  ] 
CAN  BE  EXPRESSED  AS  A  SIMPLE  LINEAR  COMBINATION  OF  THE 
CONTINUOUS  TWO-DIMENSIONAL  INVERSE  FOURIER  TRANSFORMS 
OF  THE  SPECTRUMS 


(  4  ) 


AND 


(  5  ) 


WHERE  [<l>^y],  [4>^Jy]  AND  ['I'yyl  CONTAIN  THE  (CROSS)-PSD 

FUNCTIONS  BETWEEN  THE  S  PREDICTED,  Y  OBSERVED  AND  U  TIE 
POINT  QUANTITIES. 


NOTES 

(1)  EQUATION  (3),  VIZ, 


oo  oo 

W(XQ,Yp,0)  =1  [[B(XQ-X,YQ-V,0)]-V(X,Y,H)dXdy 

— OO  — oo 

.  (  3  ) 

IS  OUR  INTEGRAL  PREDICTOR.  EACH  INDIVIDUAL  KERNEL  ELEMENT 
IN  [  B  ]  CAN  BE  THOUGHT  OF  AS  A  "WEIGHT"  SINCE  IT  IS  A 
FUNCTION  ONLY  OF  THE  DISTANCE  BETWEEN  THE  PREDICTION 
POINT  (Xo,Yo-0'  VARYING  OBSERVATIION  POINTS 

(X.Y,H=C)  HENCE,  ELEMENTS  IN  [  B  )  ARE  CIRCULARLY  SYMMETRIC 
WRT  THE  ORIGIN  AND  MONOTONICALLY  DECREASE  IN  MAGNITUDE 
AS  YOU  MOVE  AWAY  FROM  THE  ORIGIN, 

(2)  EQUATION  (  3  )  WILL  BE  NUMERICALLY  EVALUATED,  BASED  ON 
THE  FINITE  LENGTH  AND  DISCRETE  DATA  SPACING  OF  THE 
OBSERVATIONAL  GRADIENT  SURVEY  AREA. 

(3)  THE  (CROSS)-PSO  MATRICES  IN  EQUATIONS  WHS),  VIZ., 


) 


-1 

(  <1  )  * 

AND 

-1 

”  [*^yy  (^X '^y  ^  ^ ^ *^VV ^  ^ 

(  5  ) 

WILL  BE  EXPRESSED  IN  TERMS  OF 

1 )  THE  FRFQUFNCV  DOMAIN  TRANSFER  FUNCTIONS  WHICH  LINK 
THE  DISTURBING  POTENTIAL  ON  THE  GROUND,  Tq  ,  TO  THE  GROUND  ^ 

COMPONENTS  IN  M  AND  U  AND  TO  THE  AIRBORNE  GRADIENTS  IN  )£, 

2)  THE  PSD  FUNCTION  OF  Tq,  <I>^q  T0*“X'"y’’ 

ASSUME  TO  BE  ISOTROPIC,  I  E,  <1>jq  =  <I>fQ  ^q(«),  AND  IS 

BASED  ON  THE  DEFINING  PARAMETERS  OF  THE  LOCAL  GRAVITY 
MODEL  USED  AND 

3)  THE  PSD  OF  THE  WHITE  NOISE  IN  THE  GRADIENT  OBSERVATIONS 
GIVEN  BY 

T)  =  t4E2/(9/25KM2)]  =  I  nXI0'"s'‘’/(CY/M)2 


THE  NINE  TRANSFER  FUNCTIONS  THAT  ARE  APPLICABLE  TO  THIS 
STUDY  ARE: 


QUANTITY  (AT  HEIGHT  h):  F.D  TRANSFER  FUNCTION  I  INI^INO 

THE  QUANTITY  TO  Tp 


Tx 

io) 

T 

Y 

Tz 

^XX 

^XY 

Txz 

-icxi) 

Tyy 

>• 

1— 

-ioxOye"^^ 

Tzz 

WHERE  «  =  (Wj^2  +  ^,^2)I/2 

SUPPOSE  WE  WANTED  TO  PREDICT  W^q.^q)  =  (TxJyJz)  AT  EACH 

OF  THE  GGSS  GROUND  GRID  POINTS,  USING  ALL  6  GRADIENT 
OBSERVATIONS  AT  EACH  OF  THE  AIRBORNE  GRID  POINTS,  I.E, 

Y(.X,Y,H)  =  (I^XX'^XY'^XZ'^YY'^YZ'^ZZ^  tp^(c«)^,C0Y)]  WILL 

BE  THE  3  BY  6  MATRIX  HAVING  THE  STRUCTURE 


(X,XX) 

(X,XY) 

(X,XZ) 

(X,YY) 

(X,YZ) 

(X,ZZ) 

-iWw^WyN 

-«Ww^N 

-iW^y^N 

-WWj^yN 

iw^WwN 

D 

D 

0 

D 

D 

D 

(Y,XX) 

(Y,XY) 

(Y,X2) 

(Y,YY) 

(Y,YZ) 

(Y,ZZ) 

yN  -iCOj^y^N 

-tXaxWyN  -iCl)y^N 

-OXOy^N 

iw^WyN 

D 

D 

D 

D 

D 

0 

(Z,XX) 

(Z,XY) 

(Z,XZ) 

(Z.YY)  - 

(Z.YZ) 

(Z.Z2) 

(ixa^^^N 

(iKO^yN  - 

WOy^N 

-io^WyN 

-«^N 

D 

D 

D 

D 

D 

D 

WHERE  N  =  <I>^q  -^qCw)  (6) 


AND  D  =  -q  +  -  (o^^Wy^).  (7) 


(X,XX)  (X,XY) 


(X,XX) 

-iW^^N 

(X,XY) 

(X,XZ) 

-wWw^N 

(X,YY) 

-toj^y^N  • 

(X,YZ) 

-WWj^yN 

(X,ZZ) 

iw^WwN 

D 

D 

D 

D 

D  _ 

D 

(Y,XX) 

(Y,XY) 

(Y,XZ) 

(Y,YY) 

(Y,YZ) 

(Y,ZZ) 

N 

N  -oxa^yN  -i(i)y^N 

-WWy^N 

D 

D 

0 

D 

D 

D 

(Z,XX) 

(Z.XY) 

(Z,XZ) 

(Z.YY) 

(Z,YZ) 

(Z,ZZ) 

(iXO^yN 

(dOy^N 

-io^cOyN 

-W^N 

D 

D 

D 

D 

D 

D 

WHERE  N  =  4>yQ  jqCw)  (6) 

AND  D  =  T]  +  <|)jQ  •j.q(w)  (7) 

NOTES 

( 1 )  ONLY  7  OF  THE  1 6  SPECTRUMS  HAVE  TO  BE  SUBJECTED  TO  A  2D 
IFT  PROCESS  (THE  REST  ARE  DIRECTLY  ATTAINABLE  FROM  THE  7). 

(2)  DUE  TO  THE  MAKE-UP  OF  THE  DENOMINATOR  ,  D,  NONE  OF  THE  1 8 
SPECTRUMS  ARE  ISOTROPIC  WHICH  MEANS  THE  2D  IFT  PROCESS  CAN 
NOT  BE  SIMPLIFIED  INTO  A  I D  HANKEL  PROCESS. 

(3)  THE  STRUCTURE  OF  THE  SPECTRAL  MATRIX  [^.(Wx.Wy)]  WILL 
BE  A  SUBSET  OF  THE  ABOVE  [{S|]. 


C  JEKELI  (1966)  SHOWS  THAT  IF 

y  =  (■^XZ''^YZ''^ZZ^  AIRBORNE  SURVEY  POINT  AND 


W  =  T2  AT  EACH  GROUND  SURVEY  POINT 

THEN  EACH  ELEMENT  IN  THE  RESULTING  1  BY  3  SPECTRAL  MATRIX 
[p  j  ]  WILL  HAVE  THE  ISOTROPIC  DENOMINATOR  OF  - 

D  =  T]  +  ^q(O))  (8) 

AND  THE  2D  IFT  PROCESS  ON  iPi]  AND  $2]  CAN  BE  SIMPLIFIED 

TO  A  10  HANKEL  PROCESS  IN  TERMS  OF  THE  Jq  AND  J,  BESSEL 

FUNCTIONS  OF  THE  FIRST  KIND. _ 

SIMILARLY,  IF  y  =  (Txz.Tyz.T^z)  AND  W  =  (Tx,Ty,Tz)  THEN  THE 

2D  IFT  PROCESS  ON  THE  RESULTING  3  BY  3  iPj]  AND  iPj] 

MATRICES  CM  BE  SIMPLIFIED  TO  A  1 D  HANKEL  PROCESS  IN  TERMS 
OF  THE  Jq.  Jj^  and  BESSEL  FUNCTIONS _ 

NOTES: 

(1)  THE  ALGEBRAIC  MANIPULATIONS  YIELDING  THE  ID  INVERSE 
HANKEL  TRANSFORMS  ARE  VERY  METICULOUS 

(2)  CHRIS  NUMERICALLY  EVALUATED  THE  INVERSE  HANKEL 
PROCESS  VIA  SERIES  EXPANSIONS  GIVING  A  POSSIBLE  0=1  OX. 


TO  OBTAIN  N  EQUALLY  SPACED  B(X,Y)  VALUES  ALONG  EACH 
HORIZONTAL  PROFILE  (AX  =  INCREMENT)  AND  M  EQUALLY  SPACED 

B(X,Y)  VALUES  ALONG  EACH  VERTICAL  PROFILE  (AY  =  INCREMENT) 
WE  CAN  RELATE  THE  I*'*’  ANGULAR  FREQUENCIES  TO  THE  INTEGER 
FREQUENCY  COUNTERS  lj(  AND  ly  VIA 

2ji1 

Wv  = -  RADIANS /METER  AND 

NAX 


2ji-  1 

=  - -  RADIANS /METER 

M  MAY 


(15) 

AND  THEN  THE  DISCRETE  INVERSE  FOURIER  TRANSFORM  CAN  BE 
APPROXIMATED  BY 


M 

2 


B(nAX  ,mAY)=- 


MNAXAY 


2  2 

-M 

*V'2  2 


P((a^  .0)^ 

-M  1  1 


(16) 


NOTES  WRT  EQUATION  (1 6)  VIZ 


M 

2 


N 

2 


B(nAX,n]AY): 


MNAXAY 

^  Y~  2  2 


2  2 


(16). 


( 1 )  CLEARLY  THE  B  ESTIMATES  OF  ( 1 6)  APPROACH  THE  DESIRED 
CONTINUOUS  IFT  VALUES  AS  { M,N }  =>  ±  oo  AND  AS  { AX^Y)  =i^  0. 

(2)  THE  IMSL  "CANNED"  SUBROUTINE  FFT3D  CAN  COMPUTE  2D 
DISCRETE  INVERSE  FOURIER  TRANSFORMS  BY.  COMPUTING  SUMS 
OF  THE  FORM 

M-l  N-1  IL  ^  JP 

X(I+1,J  +  1)=:  22  A(L+1,P+1)©  '  N  M  ^  (17) 

?=0  L=0 


(3)  THE  “2  SIDED"  SUMS  OF  ( 16)  CAN  BE  MADE  AMENABLE  TO  THE 
IMSL  “1  SIDED"  SUMS  BY  PROPER  SHIFTING  OF  THE  SUMMATION 
OPERATORS 


(4)  FFT3D  REQUIRES  THE  USE  OF  THE  20  MATRIX  [A]  IN  ( 1 7) 
WHICH  CAUSES  STORAGE  PROBLEMS  EVEN  IF  VIRTUAL  MEMORY  IS 
INVOKED.  THE  DOUBLE  SUM  OF  ( 1 7)  CAN  BE  WRITTEN  AS 


X(I+1,J+1) 

WHERE 
G(L+1,J+1) 

(5)  ONE  CAN  A)  TRANSFORM  ALL  OF  THE  VERTICAL  PROFILES 
VIA  (19),  B)  ASSIGN  INTEGER  TAG  NUMBERS  ONLY  TO  THOSE 
RETURNING  TRANSFORMED  COMPLEX  NUMBERS  WHICH 
CORRESPOND  TO  THE  HORIZONTAL  PROFILES  OF  THE  GGSS 
SURVEY  GRID,  C)  SORT  THE  TAGGED  NUMBERS,  AND  THEN 

D)  TRANSFORM  THE  HORIZONTAL  PROFILES  VIA  ( 1 8). 

(6)  DUE  TO  THE  RADIAL  SYMMETRY  OF  THE  DESIRED  WEIGHT 
MATRIX  (  B  ],  ONLY  ONE  GGSS  QUADRANT  OF  [B]  VALUES  ARE 
NEEDED 

(7)  WITH  THE  EXCEPTION  OF  THE  SENSITIVE  B^  2Z  KERNEL,  THE 
DOMINANT  KERNEL  VALUES  NEAR  THE  ORIGIN  BEGIN  TO  CONVERGE 
WITH  “EFFICIENT"  CHOICES  OF  N,M,AX  AND  AY. 


SCENARIO; 

( 1 )  y  CONTAINS  (Tx2 JyZ'^ZZ^  SIMULATED  GRADIENTS  AT  ALL 
NODES  OF  THE  GGSS  AIRBORNE  GRID  AT  ALTITUDE  OF  H  =  BOOM. 

(2)  a  CONTAINS  PREDICTED  T2  COMPONENTS  ALONG  THE  5  N-S 
TRACKS  OF  X  =  - 1 0,-5, 0,5, 1 0  KMS ,  SPACED  EVERY  5  KMS,  FOR  ALL 
Y  £  (- 1 00  KM,  1 00  KM)  (  205  =  5(4 1 )  TOTAL  PREDICTIONS  ) 

(3)  U  CONTAINS  SIMULATED  T2  GROUND  TIE  POINT  VALUES  AS 
FOLLOWS: 


CASE  I:  2  TIE  POINTS  AT  (0,- 1 00)  AND  (0, 1 00) 

CASE  II:  3  TIE  POINTS  AT  (0,- 1 00),(0,0)  AND  (0, 1 00) 

CASE  III:  4  TIEPOINTS  AT  (-100,-1 00),(-100,1 00 ),(1 00,-1 00) 
AND  (100,100) 


MAX  ABS  ERROR 

- n 

MEAN  ERROR 

RMS  ERROR 

HANKEL  20  DISCR. 

HANKEL  2D  DISCR. 

HANKEL  2D  DISCR. 

ZAO  3.18 

0.98 

0.89 

1.39 

1.33 

2.54  2.26 

-065 

“0.47 

1.17 

0.91  ^ 

3.71  351 

113 

0.99 

1  49 

1  36 

( 1 )  CHRIS'  METHOD  OF  REDUCING  THE  AIRBORNE  GRADIENTS  IS 
QUITE  CAPABLE  OF  HANDLING  GGSS-SIZED  SURVEYS. 


(2)  THE  2D  DISCRETE  INVERSE  FOURIER  TRANSFORM  APPROACH  OF 
EVALUATING  THE  KERNELS  ALLOWS  ALL  6  GRADIENT 
OBSERVATIONS  TO  PLAY  A  ROLE  IN  THE  REDUCT  ION  PROCESS, 

(3)  ASSUMING  A  LOCAL  PSD  MODEL  IS  AVAILABLE,  THE 

PRE-DATA  REDUCTION  TASK  OF  EVALUATING  THE  KERNELS  CAN 
EASILY  BE  DONE . 

(4)  NEITHER  THE  PREDICTION  POINTS,  OBSERVATION  POINTS  NOR 
TIE  POINTS  NEED  TO  BE  REGULARLY  GRIDDED.  - 

INCONVENIENCES: 

( I )  FOR  N  =  M  =  5,056  AND  AX  =  AY  =  500  METERS,  EACH  2D 
DISCRETE  IFT  PROCESS  WOULD  TAKE  ABOUT  2  HOURS  OF  CPU  TIME 
ON  CDC/CYBER.  (THUS  FOR  EACH  SURVEY  AREA,  THE  7  NEEDED 
IFTs  WOULD  TAKE  AROUND  1 4  HOURS.) 


(2)T0  RIGOROUSLY  DETERMINE  ERROR  ESTIMATES  OF  THE 
PREDICTED  SURFACE  DISTURBANCE  COMPONENTS  REQUIRES 


[E^l  =  [Swl-2[[AJ  [A^]]. 


^  [[AJ  [A^l). 


[C^l 

tSw' 


[^YU^ 


WHERE  [A,]  =  [BlAXAY 

AND  (A2I  =  [C^]  [Cyj  +  DyJ' '  =  [P]  OF  SECTION  I. 


cm 


( I )  DUE  TO  THE  FINITE  LENGTH  OF  THE  GRADIENT  OBSERVATION 
GRID,  PREDICTIONS  NEAR  THE  PERIMETER  OF  THE  SURVEY  AREA 
WILL  BE  LESS  ACCURATE.  (CORNER  TIE  POINTS  CAN  HELP), 


(2)  A  ‘t'lQjQ  LOCAL  PSD  MODEL  MUST  BE  DEVELOPED  FOR  EACH 
SURVEY  AREA. 

(3)  ALL  AIRBORNE  OBSERVATIONS  MUST  BE  MADE  AT  THE  SAME 
HEIGHT, 


STAGE  II  SLMULATION  RESULTS  USING 
THE  NSWC  SYNTHETIC  GRAVITY  nElD 


by 

Dr.  W.  John  Hutcheson 
Bell  Aerospace  Textron 
P.O.  Box  One 
Buffalo,  NY  14240 

ABSTRACT 

The  GGSS  data  reduction  can  naturally  be  broken  down  into  two  stages. 

Stage  I,  characterized  as  being  high  frequency  and  temporal,  consists  of 
deterministic  compensations,  demodulation  and  associated  filtering.  Stage  11 
processing  Is  spatial  and  therefore  two  dimensional  in  nature  and  consists  of 
synchronous  sampling  of  the  gradients  passed  from  the  Stage  I  software,  griddlng, 
terrain  corrections,  integration,  track-crossing  adjustments,  astrogeodetic  tie 
point  adjustment,  downward  continuation  and  two  dimensional  smoothing. 

This  paper  contains  an  overview  of  the  Stage  II  algorithms  and  a  brief 
description  of  the  salient  operations  involved.  The  main  results  presented  here 
are  from  a  simulation  study  where  the  NSWC  synthetic  field  was  used  to  drive 
the  Stage  II  software.  The  effects  of  algorithm  error,  gradiometer  noise  and 
different  tie  point  configurations  are  demonstrated. 
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STAGE  I  (SEQUENTIAL)  STAGE  II  (BATCH) 


Airborne  Gradiometer  Survey 
Stage  II  Data  Reduction 


TfXTRON 


Suitability  Of  Represenatation 
Of  Integrated  Noise  Realization  By  Polynominal 

INVESTIGATION  METHOD:  SIMULATE  GGI  NOISE  ALONG  315  Km  TRACK,  PASS  THROUGH 

LSC  OPTIMAL  INTEGRATOR  AND  FIT  A  RAMP,  QUADRATIC  AND  A  CUBIC. 
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OPTIMAL  LSC  INTEGRATION  KERNEL  FOR  CENTER  POINT. 
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THE  GRADIENTS  AT  ALTITUDE  LYING  WITHIN  A  NEIGHBORHOOD  OF  THE  TRUTH  POINT  ARE  USED 
TO  TRANSFER  THE  NEAREST  MAP  VALUE  OF  T.  DOWN  TO  THE  TRUTH  POINT. 


Central  Area  Of  NSWC  Synthetic  Field 
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TRUTH/ESTIMATE  COMPARISON 


CENTER  EAST/WEST  TRACK  TX 
TRUTH  AND  ESTIMATE 


Overall  Error  Standard  Deviation 
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CENTRAL  EAST/WEST  TRACK  TX 
TRUTH  AND  ESTIMATE 


Main  Results 
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TITLE  OF  PAPER:  Stage  II  Simulation  Results  Using  the 
MSWC  Synthetic  Gravity  Field 

SPEAKER;  A1  Jlrcltano 

QUESTIONS  AND  COMMENTS: 

1.  Question ;  Chris  Jekeli 

What  was  the  white  noise  of  the  gradiometer? 

Response: 

50  e2/H2 

2.  Question :  Alan  Rufty 

Were  all  tracks  of  predicted  answers  coincident  with  the  given  data  tracks? 
Response: 

Yes,  but  the  answers  degrade  minimally  as  one  goes  away  from  the  raw  data 
tracks 

3.  Question :  Richard  Rapp 

What  was  the  accuracy  of  the  tie  points? 

Response ; 

Approximately  0.1  arcsec  for  deflections  and  0.1  mgal  for  disturbances. 

4.  Question ;  John  Brozena 

Are  the  tie  point  data  used  to  constrain  the  least  squares  track 
adjustments? 

Response : 


Yes. 
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SPHERICAL  HARMONIC  EXPANSIONS 

by 
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Ohio  State  University 
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Golumbus,  Ohio  43210-1247 

ABSTRACT 

Spherical  harmonic  expansions  can  be  used  to  describe  the  earth's 
gravitational  field.  The  resolution  of  these  fields  depends  on  the  highest 
degree  in  the  field.  In  the  past  year  several  fields  to  degree  180  (0SU81), 
200  (GPM2),  250  (OSU86C/D)  and  360  (0SU86E/F)  have  become  available.  Such 

fields  are  needed  to  compute  geoid  undulations,  deflections  of  the  vertical, 
etc.,  or  to  provide  a  reference  field  for  reductions  of  local  data,  such  as 
altimeter  or  gradiometer  data.  This  presentation  will  consider  the  gradient 
information  in  these  new  fields  and  will  compare  solution  differences  with 
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Abstract 


Spherical  harmonic  expansions  can  be  used  to  describe  the  earth’s 
gravitational  field.  The  resolution  of  these  fields  depends  on  the  highest 
degree  in  the  field.  In  the  past  year  several  fields  to  degree  180  (0SU81), 
200  (GPM2},  250  (OSU86C/0)  and  360  (0SU86E/F)  have  become  available.  Such 
fields  are  needed  to  compute  geoid  undulations,  deflections  of  the  vertical, 
etc.,  or  to  provide  a  reference  field  for  reductions  of  local  data,  such 
as  altimeter  or  gradiometer  data.  This  presentation  will  consider  the 
gradient  information  in  these  new  fields  and  will  compare  solution 
differences  with  formal  accuracy  estimates  to  assess  the  accuracy  of  these 
new  fields. 


High  Degree  Fields 


Recent  Developments 


•  OSU81 

•  GPM2  - 1985 

•  OSU86C/D 


•  OSU86E/F 


•  June  1986  1  x  1  Terrestrial  Anomalies 

•  1985  1  X  1  Altimeter  Derived  Anomalies 

•  GEML2’  Potential  Coefficients 

•  OSU86D  Uses  Geophysically  Predicted 
Anomalies 

•  OSU86C  Excludes  Geophysically 
Predicted  Anomalies 

•  Least  Squares  Combination  Followed  By 
Rigerous  Optimal  Estimation  To  n:250 


New  High  Degree 

Fields 

OSU86E,  OSU86F 

•  August  1986  30’x30’  Terrestrial  Anomalies 

•  1985  30x30’  Altimeter  Derived  Anomalies 

•  Solutions  Made  By  Forcing  Mean  of 
30’x30’  Values  to  Agree  with  Adjusted  1 
Values 

•  OSU86E  -  No  Geophysical  Anomalies 

•  OSU86F  -  Includes  Geophysical 
Anomalies 

•  Coefficient  to  n:360  by  Quadratures 
(HARMIN) 


Solution  Comparisons 


•  Accuracy  Estimates 

•  Anomalies 

•  Gravity  Disturbances 

•  Deflections  of  the  Vertical 


Gradients 


RMS  Values  implied 
By  OSU86F  Field 


To  180  To  360 

Anomaly  (mgal) 
Disturbance  (mgal) 
Deflection  (secs) 
Undulation  (m) 
Gradient  (Tzz)O 
Gradient  (Txx)(E) 

♦  24.  6  t26.9 
i  30.  3  t32.3 
t  6.  0  i  6.4 
t  30 .  4  i304 
t  3.  4  i5.5 
t  1.7  -^2.8 

*  on  the  surface  of  a  sphere  of  radius 
6371  km 


Comparison  of  RMS  Values 
Implied  By  OSU86F  and 
GPM2  to  Degree  180 


OSU86F 

GPM2 

Anomaly  (mgal) 

t  24. 

6 

t26.8 

Disturbance  (mgal) 

t  30. 

3 

t32.2 

Deflection  (secs) 

+  6. 

0 

+  6.4 

Undulation  (m) 

t  30. 

4 

♦305 

Gradient  (Tzz)(E) 

t  3. 

4 

♦4.1 

Gradient  (Txx)(E) 

t  1  . 

7 

♦21 

*  on  the  surface  of  a  sphere  of  radius 
6371  km 


Applications 
in  Gradiometry 


Provide  Long  Wavelength  Information 


•  Definition  of  Long  Wavelength  (>500  km) 

•  Corresponding  Degree  About  80 

•  Gravity  Disturbances  of  Prime  Interest? 

•  Accuracy  of  OSU86D 

•  Comparison  of  OSU86F  and  GPM2 


Anomaly  Degree  Variances  and  Their  Accuracy 
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Gravity  Disturbance  Accuracy:  OSU86D  solution 
Gravity  Disturbance  Difference:  OSU86F  -  GPM2 
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Gravity  Disturbance  Accuracy:  OSU86D  solution  ^ 
Gravity  Disturbance  Difference;  OSU86F  -  GPM2 
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Deflection  Accuracy:  OSU86D  solutions 
Deflection  Difference:  OSU86F  -  GPM2 


30.00  60.00  90.00  120.00 

HARMONIC  DEGREE 


150.00 


180.00 


•  Improved  High  Degree  Fields  Exist 

•  RMS  Disturbance  Accuracy  is  4.3  mgal 

•  RMS  Disturbance  Difference  in  the  U.S. 
Between  OSU86  and  GPM2  is  2.5  mgal 

•  RMS  Global  Disturbance  Difference  is 
4.8  mgal 

•  RMS  Total  Deflection  Accuracy  is  0.9  secs 

•  RMS  Total  Deflection  Difference  is  1.0  secs 

•  All  above  values  for  degrees  2  to  80 


Gravity  Disturbance  Differences  (OSU86F-GPM2) 

1  mgal  contour  interval 
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TITLE  OF  PAPER: 


Gradient  Information  in  New  High  Degree  Spherical  Harmonic 
Expansions 


SPEAKER:  Richard  H.  Rapp 
QUESTIONS  AND  COMMENTS: 


1.  Question:  Charles  F.  Martin 

To  what  extent  is  the  observational  data  powerful  enough  to  support 
harmonic  expansion  up  to  degree  and  order  360“? 

Response : 


It  depends  on  the  quality  of  data  over  the  areas  of  interest,  i.e.,  U.S., 
Central  Europe,  Marianas  Trench.  Input  data  of  high  quality  and  harmonic 
field  above  180“  provides  significant  information. 

2.  Question:  A1  Jlrcltano 

Is  accuracy  of  gravity  data  better  in  ocean  areas  or  land  areas? 

Response : 


Generally  better  in  land  areas.  In  the  US  and  Europe,  the  accuracy  is  about 
2-3  mgal  in  1“  squares;  in  ocean  areas  it  is  about  6-7  mgal  in  1“  squares. 

3.  Question:  Warren  Heller 

What  do  you  see  as  the  primary  error  sources  driving  the  approximately  4 
mgal  of  error  in  harmonics  through  degree  80? 

Response: 

Primarily  surface  data  quality  and  data  omission.  (Some  discussion  of 
advantages  of  considering  local  areas  where  data  is  good;  also  discussion  of 
correlation  between  satellite  and  terrestrial  gravity  measurement  errors). 

4.  Question :  Jim  Lowrey 

Are  there  any  plans  to  extend  the  model  out  beyond  360“? 

Response: 


Curently  there  are  none;  however,  an  extension  out  to  720 
would  be  possible  although  the  need  for  this  is  questionable. 


